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FOREWORD 

This  document  presents  the  results  of  a  study  of  improvements  which 
could  be  made  in  aerodynamic  heating  programs  to  make  them  better  able  to 
handle  complex  high  speed  vehicles.  The  study  was  conducted  by  the  High 
Speed  Aero  Performance  Branch  (FXG) ,  Aeromechanics  Division,  Air  Force 
Flight  Dynamics  Laboratory,  Wright-Patterson  AFB,  Ohio.  The  work 
concludes  an  in-house  program  and  was  performed  under  task  240407 
"Aeroperformance  and  Aeroheating  Technology",  work  unit  24040714, 

"Development  of  Numerical  Techniques  for  Predicting  Aerodynamic  Heating 
to  Flight  Vehicles"  and  covers  work  performed  from  May  1975  to  September  1977. 


Ac  a  ess  ion  For 


X 


NlIS  GBlA&I 
DtDC  TAB 
Unannounced 
Justification 


By 


Distribution/ 


Availability  Codes 


Dist. 


ft 


Avail  and/or 
special 


D  D  C 


MZL1? 

TT\ 

W  DEC  12  1979 

I,,. 

JiiEEEinnp 

—> 

D 


iii 


'  S' v 


AFFDL- TR- 79-3001 


TABLE  OF  CONTENTS 


SECTION  PAGE 

I  INTRODUCTION  1 

II  MARK  IV  SUPERSONIC-HYPERSONIC  ARBITRARY  4 

BODY  PROGRAM  STUDIES 

1.  Program  Description  4 

2.  Blunt  Cone  Studies  with  the  Mark  IV  HABP  5 

3.  X-24C-10D  Studies  with  the  Mark  IV  HABP  6 

4.  Summary  of  the  Mark  IV  HABP  Studies  8 

III  STUDIES  WITH  THE  DEJARNETTE  PROGRAM  10 

1.  Program  Description  10 

2.  Program  Limitations  and  Sources  of  Elvov  15 

3.  Sample  Test  Cases  18 

IV  THREE-DIMENSIONAL  SHOCK-CAPTURING  TECHNIQUE  27 

1.  Background  27 

2.  Cone  Study  29 

3.  Ogive-Cylinder  Study  31 

4.  Lifting  Body  Study  33 

V  HEAT  TRANSFER  CALCULATIONS  USING  SCT  PRESSURES  39 

AND  DEJARNETTE  PROGRAM 

VI  CONCLUSIONS  42 

REFERENCES  44 

APPENDIX  A:  BODY  PROGRAM  FOR  OGIVE-CYLINDER  49 

APPENDIX  B:  BODY  PROGRAM  FOR  A  GENERAL  BODY  51 


v 


nrjrminwi  i 


' 1*  V; ..  YV. •. 


AFFDL-TR-79-3001 


LIST  OF  TABLES 


TABLE  PAGE 

I  Data  for  Six  Sample  Test  Cases  58 

II  Comparison  of  Desired  and  Smoothed  Phi  Angles  59 

III  Location  of  Circumferential  Planes  for  X-2AC-10D  59 

Cases 


vi 


AFFDL-TR-79-300'1 


LIST  OF  ILLUSTRATIONS 

FIGURE  TITLE  PAGE 

1  X-24C-10D  Geometry  60 

2  Comparison  of  Numerical  and  Experimental  62 

Stanton  Numbers  for  the  X-24C-10D  Side 

Panel  at  Model  Station  20.52 

3  Location  of  X-24C-10D  Streamlines  (HABP)  64 

4  Plot  of  Streamline  Angle  Vs  Angle  of  Attack  65 

5  Orthogonal  Coordinate  System  and  Body  66 

Orientation  for  the  DeJarnette  Program 

6  25°  Cone  Pressure  Distribution  67 

7  25°  Cone  Heat-Transfer  Distribution  68 

8  0. 7:1.0  Elliptical  Cone  Description  69 

9  0.7: 1.0  Elliptical  Cone  Windward  Streamline  70 

Heat-Transfer  Distribution 

10  0. 7:1.0  Elliptical  Cone  Circumferential  71 

Heat-Transfer  Distribution 

11  80°  Delta  Wing  Windward  Streamline  Heat-  74 

Transfer  Distribution 

12  80°  Delta  Wing  Spanwise  Heat-Transfer  Distribution  75 

13  10°  Cone  Windward  Streamline  Heat-Transfer  76 

Distribution 

14  10°  Cone  Spanwise  Heat-Transfer  Distribution  77 

15  Ogive-Cylinder  Geometry  79 

16  Surface  Pressures  on  Ogive-Cylinder  at  Model  80 

Station  41.7 

17  Pitot  Pressure  Surveys  for  the  Ogive-Cylinder  at  81 

Model  Station  41.7 

18  Surface  Pressures  on  the  X24C-10D  at  2°  Angle  84 

of  Attack 

19  Surface  Pressures  on  the  X24C-10D  at  6°  Angle  of  87 

Attack 

vii 


AFFDL-TR-79-3001 


LIST  OF  ILLUSTRATIONS  (Coat'd) 

Figure  Title  Page 

20  X24C-10D  Inviscid  Shock  Location  (Longitudinal  90 

Plane) 

21  X24C-10D  Inviscid  Shock  Location  (Tangential  Plane)  91 


22  Longitudinal  Slope  of  the  Shock  Wave  at  $  =  0°  93 

23  Angular  Dependence  of  the  Longitudinal  Slope  of  94 

the  Shock  at  Three  Body  Stations 

24  Tangential  Slope  of  the  Bow  Shock  95 

25  Comparison  of  Numerical  and  Experimental  97 

Stanton  Numbers  on  the  Ogive-Cylinder 

26  Comparison  of  Numerical  and  Experimental  Heat  101 

Transfer  to  the  Ogive-Cylinder 

27  Comparison  of  Numerical  and  Experimental  102 


Pressures  and  Stanton  Numbers  on  the 
Ogive-Cylinder  as  a  Function  of  Angle  of 
Attack 


A-l  Construction  Used  in  the  Calculation  of  the  103 

Ogive-Cylinder  Geometry 

B-l  Typical  Cross  Section  of  the  X24C-10D  104 

B-2  Construction  of  the  Body  Radius  for  a  105 

Circular  Bottom 

B-3  Construction  of  the  Body  Radius  for  a  General  106 

Circular  Segment 

B-4  Construction  Used  for  a  Straight  Bottom  Segment  107 

B-5  Construction  Used  for  a  General  Straight  Segment  108 


viii 


A*  .  ■  ■  . 

■.W4-*':  -- 


AFFDL-TR-79-3001 


LIST  OF  SYMBOLS 

See  Eq  B-l 
See  Eq  B-l 
See  Eq  B-l 

See  Eq  B-24  or  Figure  B-5 

See  Eq  B-l 

See  Figure  B-2 

See  Eq  B-l 

See  Eq  B-l 

Body  Radii  used  in  four-point  numerical 
differentiation  formula.  See  Eqs  B-31, 
32,  and  33 

Local  coefficient  of  skin  friction 
Radius  of  the  bow  shockwave 
See  Eq  B-l 
See  Eq  B-l 


Cy  See  Eq  B-l 

Dtemp  See  Eq  A— 5 

f  In  the  DeJarnette  analysis,  body  radius  measured  normal 

to  the  longitudinal  body  axis.  See  Figure  5 

h  Either  body  scale  factor  or  local  heat-transfer 

coefficient  in  Btu/f t2-sec-°R.  In  Eqs  B-31,  32  and 
33  the  interval  between  and  B2. 

ho2  Stagnation  point  heat-transfer  coefficient  (Btu/f t2-sec-°R) 

H  Enthalpy 

Hw  Wall  enthalpy 


ix 


AFFDL-TR- 79-3001 


Ho 

K 

L 

Log 

M 

P 

n 

r 

R 

H 

^Base 

Hey 

Hn 

Hog 

Rtemp 


Ho 

R1  ‘ 


Reoo 

S 


sog 

Stgo 


Total  enthalpy 

Curvature  of  the  boundary  layer  edge  streamline  in 
the  plane  tangent  to  the  body  surface.  Also  used 
as  a  dummy  variable  in  Appendix  B 

Body  length.  See  Table  1 

Length  of  blunted  ogive.  See  Figure  A-l 

Mach  number 

Pressure 

Local  aerodynamic  heat-transfer  rate 
Length  in  the  radial  direction 
Radius 

Local  body  radius.  See  Figure  B-l 
Base  radius.  See  Page  29 
Radius  of  cylinder.  See  Figure  A-l 
Nose  radius.  See  Page  29  U  Table  I 
Radius  of  ogive.  See  Figure  A-l 
Dummy  variable  defined  by  Eq  A-5 

Radii  of  circular  arcs.  See  Figure  B-l 


Freestream  unit  Reynolds  number  Cl/ ft) 

In  the  DeJarnette  analysis,  S  is  the  distance  along  a 
streamline  measured  from  the  stagnation  point 

Length  of  ogive.  See  Figure  A-l 

Stanton  number  based  on  freestream  conditions 
and  freestream  total  temperature 


X  Temperature 

u  In  the  DeJarnette  analysis  the  velocity  in  the 

streamline  direction;  in  the  NSWC  analysis,  the 
tangential  velocity 


x 


AFFDL-TR-79-30G1 


v 

V 

voo 

W 

Wf 

X 


xa 

Xn 


xtre 


xvn 


Xo 

*1 


Y 


Ya 


Velocity  in  the  radial  direction 
Frees tream  velocity 
Velocity  In  the  axial  direction 
Weighing  factor.  See  Eq  (4) 

In  the  DeJarnette  analysis,  distance  along  the 
longitudinal  axis  of  the  body  from  the  nose- tip 
See  Figure  1  for  coordinate  system  description 

See  Figure  B-l 

See  Eq  B-l 

In  the  DeJarnette  analysis,  the  distance  along  the 
longitudinal  axis  of  the  body  between  the  nose- tip 
and  the  end  of  boundary  layer  transition 

In  the  DeJarnette  analysis,  the  distance  along  the 
longitudinal  axis  of  the  body  between  the  nose- tip 
and  the  beginning  of  boundary  layer  transition 

In  the  DeJarnette  analysis,  the  distance  along  the 
longitudinal  axis  of  the  body  measured  from  the  vertical 
origin 


Coordinates  used  in  the  construction  of  the 
X24C-10D  cross  section.  See  Figure  B-l 


In  the  DeJarnette  analysis,  the  distance  normal  to 
the  longitudinal  axis  of  the  body  in  the  X-Y  plane. 
See  Figure  1  for  coordinate  system  description 

See  Figure  B-l 

See  Eq  B-l 


Coordinates  used  iu  the  construction  of  th£ 
X-24C-10D  cross  section.  See  Figure  B-l 


In  the  DeJarnette  analysis,  the  distance  normal  to  the 
longitudinal  axis  of  the  body  in  the  X-Z  plane.  See  Figure 
1  for  coordinate  system  description.  In  the  NSWC  analysis 
the  longitudinal  coordinate  in  the  cylindrical  coordinate 
system 


xi 


AFFDL-TR-79-3001 


Zmax  In  the  DeJarnette  analysis  the  maximum  value  of  Z  for 

a  particular  cross  section 

a  Angle  of  attack.  See  Figure  5.  Also  in  the  NSWC  geometry 

analysis  a  body  scope  at  sphere-ogive  function.  See  Figure  A-l 

y  Ratio  of  specific  heats.  In  the  NSWC  analysis  a  body  angle. 

See  Figure  B-3 

6  Boundary  layer  displacement  thickness 

Q  Angle  used  in  body  description.  See  Figures  B-2  and  B-4 

0nl  Boundary  layer  momentum  thickness 

y  In  the  DeJarnette  analysis  the  viscosity.  In  the  NSWC 

analysis  a  substitution  used  in  the  body  description. 

See  Eq  B-14. 


p  Density 

a  A  substitution  used  in  the  body  description 

See  Eq  B-15.  Also  slope  of  a  straight  line  segment  in 
X-24C-10D  cross  section.  See  Figure  B-l 

4>  Circumferential  body  angle.  See  Figures  1  and  B-l 

4>m  Angle  used  in  body  description.  See  Figure  B-3 

4>s  Angle  used  in  body  description.  See  Figure  B-5 

Angle  used  in  body  description.  See  Figures  B-2  and  B-3 

Subscripts 

AW  Adiabatic  wall  conditions 

b  Body 

cy  Cylinder 

e  Conditions  the  inviscid  flow  at  the  boundary  layer  edge 

o  Total  or  reservoir  conditions,  or  a  given  value 

og  Ogive 

02  Conditions  at  the  stagnation  point 

W  Conditions  at  the  wall 

®>  Freestream  conditions  in  the  undisturbed  flowfield 

upstream  of  the  vehicle 

xii 


•  -.V 


AFFDL-TR-79-3001 


SECTION  I 
INTRODUCTION 

The  development  of  complex  lifting  configurations  and  high  speed 
maneuvering  vehicles  has  emphasized  the  need  for  numerical  techniques 
to  predict  aerodynamic  heating  rates  as  a  function  of  the  vehicle 
trajectory.  These  numerical  programs  are  not  expected  to  eliminate 
the  requirement  for  wind  tunnel  and  flight  testing,  but  are  expected 
to  result  in  efficient  use  of  test  time  and  improve  confidence  that 
all  potential  problem  areas  on  the  vehicle  have  been  examined. 

Large-scale  numerical  programs  for  the  prediction  of  force  and 
moments  are  available  for  complex  shapes.  One  such  program,  the  Supersonic- 
Hypersonic  Arbitrary  Body  Program  (HABP)  (Reference  1) ,  also  includes 
routines  to  trace  streamlines  and  calculate  viscous  effects  including 
heat  transfer  rates  to  the  body.  An  earlier  version,  the  Mark  III 
Hypersonic  Arbitrary  Body  Program  (HABP)  (Reference  2) ,  contained  a 
simplified  skin  friction  option  to  calculate  viscous  effects.  The  Mark 
III  version  of  the  program  has  been  modified  and  used  extensively  at 
NASA  Langley  to  predict  heating  rates  on  shuttle-type  vehicles  at  many 
points  along  the  reentry  trajectory  (Reference  3).  While  the  Mark  IV 
HABP  Still  contains  the  simplified  skin  friction  option,  it  also  has  an 
option  to  treat  viscous  effects  with  an  integral  boundary  layer  method. 

Very  little  experience  has  been  documented  with  the  integral  boundary 
layer  option  in  the  Mark  IV  HABP. 

A  second  program  to  predict  the  aerodynamic  heating  to  shuttle-type 
vehicles  has  been  developed  by  DeJarnette  in  cooperation  with  NASA  Langley 
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(Reference  4) .  This  program  is  based  on  the  axir;  mmetric  analog  and  has 
been  modified  to  include  entropy  swallowing,  input  pressure  distributions 
and  a  more  general  body  geometry.  Although  basically  an  ideal  gas 
program,  it  does  offer  an  equilibrium  air  option. 

A  third  program,  named  MINIVER,  was  developed  by  the  McDonnell 
Douglas  Company  (Reference  5).  MINIVER  is  a  smaller  version  of  a  more 
complex  heating  program  and  differs  significantly  from  the  HABF  and 
DeJarnette  programs  in  that  it  is  a  point-calculation  program.  In  such  a 
program,  the  local  pressure  is  found  and  a  running  length  from  the 
stagnation  point  to  the  point  in  question  is  approximated  in  some  manner; 
a  flat  plate  analogy  is  used  in  the  present  code.  One  major  advantage 
of  this  code  is  that  the  heat  transfer  rates  can  be  obtained  at  a  set  of 
points  which  are  input  by  the  user.  This  advantage  is  off-set  by  the 
need  to  know  the  origin  and  length  of  the  streamline  croasing  that  point. 

The  HABP  was  used  to  calculate  heat  transfer  rates  on  the  X-24C-10D 
configuration;  the  DeJarnette  code  was  used  with  the  ogive-cylinder.  The 
MINIVER  code  was  not  used,  except  to  demonstrate  proper  operation  on  the 
local  computer  system,  due  to  the  lack  of  streamline  information  for 
complex  vehicles.  The  MINIVER  program  may  prove  to  be  very  useful  in  the 
future  for  use  on  wings,  fins  and  possibly  the  side  panels  of  an  X-24C 
type  configuration. 

During  this  study  it  became  apparent  that  an  improvement  in  program 
accuracy  could  be  obtained  with  a  better  Inviscid  surface  pressure 
calculation.  At  about  the  same  time,  the  Three-Dimensional  Shock- 
Capturing  Technique  (SCT)  based  on  work  by  Kutler  and  developed  by  Solomon 
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(Reference  6)  became  available.  This  code  proved  to  be  fast  and  reliable, 
but  lacked  an  adequate  geometry  routine  for  the  type  of  bodies  of  interest. 
New  geometry  routines  were  developed  for  the  SCT  and  one  case  was  run 
using  the  SCT  pressure  distribution  in  the  DeJarnette  heat  transfer 
program. 

During  the  period  these  numerical  methods  were  being  studied, 
experimental  investigations  of  some  configurations  of  interest  were 
also  being  conducted  at  the  Arnold  Engineering  Development  Center  by  AFFDL. 
The  experimental  results  were  intended  for  comparison  with  the  numerical 
results.  There  is  never  enough  experimental  data  to  really  satisfy  the 
comparison  requirements  of  the  numerical  program;  in  this  case,  hindsight 
indicates  more  pitot  pressure  surveys,  flow  angularity  studies  and 
better  knowledge  of  the  location  of  both  bow  and  imbedded  shocks  would 
have  increased  confidence  in  the  numerical  technique.  A  follow-on  study 
will  address  these  problems. 

Some  of  the  comparisons  are  adequately  covered  in  other  reports  and 
are  only  referenced  here.  The  purpose  of  this  report  is  not  to  detail 
all  of  the  results  developed  by  this  study,  but  to  indicate  thi  present 
capabilities  of  the  various  numerical  techniques  and  Indicate  areas  in 
which  improvements  have  been  or  can  be  made. 

The  ultimate  goal  of  this  effort  was  to  extend  the  aero  heating 
prediction  capabilities  of  HASP,  the  DeJarnette  code  and  the  inviseid 
pressure  methods  to  practical  bodies.  No  improvements  were  possible  with 
the  HABP  at  this  time,  but  significant  improvements  were  made  in  the 
DeJarnette  program  and  the  inviseid  pressure  program. 
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SECTION  II 

MARK  IV  SUPERSONIC-HYPERSONIC  ARBITRARY  BODY  PROGRAM  STUDIES 
1.  PROGRAM  DESCRIPTION 

The  Mark  IV  HABP  was  developed  by  Gentry  e_t  al.  at  the  Douglas 
Aircraft  Division  of  the  McDonnell  Douglas  Corporation  under  contract 
to  the  AF  Flight  Dynamics  Laboratory  (Reference  1).  Earlier  versions  of 
the  program  were  developed  with  support  from  Douglas  Independent  Research 
and  Development  funds  and  Air  Force  contracts  (Reference  2).  One  goal  of 
the  program  was  to  supply  design  engineers  with  a  flexible  engineering 
tool  to  provide  realistic  modeling  of  the  flow  about  complex  aerodynamic 
shapes.  While  designed  primarily  to  supply  first-order  aerodynamic  data, 
the  program  also  includes  viscous  two-dimensional  boundary  layer 
calculations.  A  wide  choice  of  empirical  pressure  techniques  is  available 
for  selection  by  the  user;  the  resultant  pressure  distribution  is  used  to 
obtain  the  inviscid  flow  field. 

The  Mark  III  HABP  contains  a  simplified  skin  friction  option.  The 
Mark  XV  HABP,  in  addition  to  this  skin  friction  option,  contains  an 
integral  boundary  layer  technique  based  on  McNally  (Reference  7)  which  was 
the  option  used  in  this  study.  To  use  this  option,  the  program  is  first 
required  to  locate  the  streamline  trajectory  along  the  body  surface. 

The  calculation,  although  carried  out  on  a  three-dimensional  body,  is 
basically  a  two-dimensional  calculation;  that  is,  cross-flow  pressure 
gradients  are  neglected.  The  integral  boundary  layer  method  uses  Cohen 
and  Reshotko's  (Reference  40)  method  £or  laminar  boundary  layer  calculations, 
Sasman  and  Cresci'a  (Reference  41)  method  for  turbulent  boundary  layers  and 
the  Schlichting-Wrich-Granville  method  for  predicting  the  transition 
point.  The  present  version  does  not  calculate  transitional  flow;  it 
switches  from  laminar  to  turbulent  at  the  transition  point. 
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The  Mark  IV  HABP  obtains  its  geometry  flexibility  by  allowing  the 
user  to  construct  the  body  as  a  family  of  flat  elements.  Surface  elements 
are  combined  to  form  a  section  of  the  vehicle,  sections  are  combined 
into  panels,  and  panels  combined  into  vehicle  components.  Regardless  of 
what  type  of  geometry  input  is  used,  eventually  a  set  of  plane, 
quadrilateral  surface  elements  are  generated  and  their  centroids  stored. 

All  surface  data  are  calculated  based  on  these  centroids.  Each  streamline 
calculation  begins  at  the  centroid  of  a  surface  element  and  continues 
to  the  centroid  of  the  last  element  of  the  panel  on  which  it  started. 

There  is  no  calculation  of  streamlines  from  one  panel  to  another. 

2.  BLUNT  CONE  STUDIES  WITH  THE  MARK  IV  HABP 

The  first  cases  attempted  with  the  Mark  IV  HABP  were  blunted  cones. 

In  addition  to  the  obvious  analytic  advantages  of  studying  a  cone  first, 
the  cone  is  a  particularly  easy  body  to  input  into  the  program.  The  two 
cones  selected  were  a  10°  cone  with  a  1-inch  nose  radius  at  Ha  ■*  7.98 
and  a  25°  cone  with  a  0.5— inch  nose  radius  at  M®  •  7.95.  Successful 
solutions  were  not  obtained  for  either  cone  at  any  angle  of  attack. 

Several  factors  contributed  to  problems  with  the  cone  solutions,  but  the 
most  important  problem  was  the  inability  of  the  program  to  obtain  a  successful 
surface  fit  of  the  pressure  distribution. 

To  understand  this  failure,  it  is  necessary  to  know  that  the  pressure 
at  each  element  centroid  is  fit  with  a  surface  spline.  While  the  cone  is 
an  easy  analytic  problem,  the  resultant  pressure  distribution  is  not  simple 
if  the  entire  cone  is  to  be  treated  as  one  panel.  Successful  fits  could 
be  obtained  by  dividing  the  cone  into  several  panels  and  fitting  the  pressure 
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on  each  panel  individually,  but,  of  course,  if  this  were  done,  the 
streamlines  could  not  be  traced  along  the  cone.  Further  efforts  with 
a  cone  were  discontinued.  In  closing  this  portion  of  the  work,  it 
should  be  noted  the  Mark  IV  HABP  was  never  intended  to  perform  this 
type  of  calculation. 

3.  X-24C-10D  STUDIES  WITH  THE  MARK  IV  HABP 

When  experimental  pressure  data  became  available  from  the  X-24C-10D 
this  configuration  was  programmed  into  the  Mark  IV  HABP.  Figure  la*  shows 
the  model  configuration;  Figure  lb  the  simplified  side  view  used  in  the 
HABP.  Initial  efforts  resulted  in  problems  similar  to  those  encountered 
with  the  cone.  The  problem  was  resolved  by  making  the  flat  bottom 
surface  on  the  vehicle  one  panel  and  the  upper  surface  a  second  panel. 

Not  all  details  were  included  in  the  numerical  model.  Because  the  surface 
pressures  in  the  HABP  are  predicted  using  empirical  methods  (i.e.  the 
surface  pressure  is  predicted  as  a  function  of  the  turning  angle  and 
freestream  raach  number),  the  bottom  panel  was  a  constant  pressure  surface 
and  did  not  even  approximate  proper  streamline  locations  or  viscous 
effects.  The  HABP  does  permit  the  input  of  surface  pressures  at  the 
centroids  of  each  panel,  but  adequate  surface  pressures  were  not  available 
If  the  surface  pressures  were  known,  empirical  methods,  for  example 
DeJarnette,  are  available  for  calculating  the  viscous  effects  on  such  a 
surface  with  the  same  accuracy  of  the  HABP. 

However,  an  extensive  series  of  computer  runs  was  made  on  the  side  of 
the  X-24C-10D.  Experimental  Stanton  Numbers  were  available  (Reference  8) 
at  station  20.52  on  the  model.  In  the  HABP,  streamline  calculations  must 

^Figures  and  tables  are  located  at  end  of  report,  pages  58  through  108. 
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start  at  the  centroid  of  an  element.  As  none  of  the  streamlines  starting 
at  centroids  near  the  nose  remained  on  the  model  side  at  this  station, 
it  was  necessary  to  start  streamlines  at  centroids  on  the  leading  edge 
of  the  model.  Different  centroids  were  selected  for  each  angle  of 
attack  calculated,  moving  further  back  with  increasing  angle.  In  general, 
good  agreement  was  obtained  with  experimental  results  as  shown  in  Figure  2. 
There  were  however,  several  problems  which  should  be  noted  by  future 
users.  The  Mark  IV  HABP  locates  the  streamlines  by  an  integration  of 
the  Euler  equations.  These  steps  may  be  quite  small  (relative  to  the 
length  of  the  panel)  particularly  in  regions  of  steep  pressure  gradients. 

As  storage  locations  are  available  for  only  100  streamline  points,  the 
user  may  choose  to  save  only  every  nfc^  point.  This  resulted  in  the  saving 
of  only  every  5t^1  point  on  many  streamlines  and  occasionally  only  every 
10th  point  on  long  streamlines  passing  through  regions  of  steep  pressure 
gradients.  This  restriction  resulted  in  warnings  from  the  integral 
boundary  layer  program  that  the  pressure  change  between  data  points 
exceeded  the  values  allowed  by  the  tables  used  in  the  program.  Although 
there  was  no  incidence  of  program  failure  due  to  this,  and  the  resultant 
heat  transfer  rates  compare  favorably  with  experimental  data,  such  warning 
can  not  be  ignored  and  all  such  solutions  must  be  treated  with  suspicion. 

Three  changes  might  be  made  to  correct  this  problem.  New  tables  should 
be  included  in  the  integral  boundary  layer  routine  to  permit  steeper 
pressure  gradients,  the  code  could  be  changed  to  permit  an  input  value  to 
change  the  points  saved  in  a  particular  region  of  the  body,  or  more  storage 
could  be  set  aside  to  store  streamline  data.  While  the  additional  storage 
would  appear  to  be  the  easiest  solution,  the  elaborate  mass  storage  file 
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built  by  the  Mark  IV  HABP  discouraged  any  attempt  at  restructuring.  While 
the  other  two  are  possible  modifications  to  the  code,  a  decision  was  made 
not  to  change  the  program  at  this  time. 

Because  the  empirically  predicted  pressure  is  constant  on  a  flat 
surface,  the  streamlines  calculated  by  the  Mark  IV  HABP  are  parallel 
straight  lines  (Figure  3)  on  the  flat  side  of  the  vehicle.  Some  curvature 
can  be  =.nen  on  the  crown.  To  determine  how  representative  these  streamlines 
really  are,  approximate  flow  angularity  was  obtained  from  oil  flow 
experiments  and  compared  with  the  numerical  results.  As  seen  in  Figure  4, 
the  results  are  acceptable.  A  single  streamline  angle  from  Powers'  method 
of  characteristics  program  (Reference  31)  is  also  shown. 

4.  SUMMARY  OF  THE  MARK  IV  HABP  STUDIES 

The  Mark  IV  HABP  was  found  to  be  a  reasonably  accurate  program  for  the 
prediction  of  viscous  effects  on  general  shapes.  If  the  program  is  to  be 
used  on  a  wide  variety  of  shapes,  tables  used  in  the  integral  boundary  layer 
routine  should  be  extended  to  handle  a  wider  range  of  pressure  gradients, 
the  number  of  streamline  points  saved  should  be  increased  or  the  number  of 
points  saved  should  be  a  function  of  pressure  gradient.  The  major  advantage 
of  this  program  is  its  very  flexible  geometry  package;  this  advantage  is 
off-set,  however,  by  the  current  restriction  on  saving  streamline  data 
and  the  inability  to  continue  streamlines  across  panel  boundaries. 

It  is  important  to  realize  the  selection  of  panels  on  the  vehicle  may 
be  governed  by  the  behavior  of  the  surface  spline  used  in  the  program. 

For  a  surface  spline  to  result  in  acceptable  fits  the  region  fit  must  be 
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normalized  in  a  manner  consistent  with  the  behavior  of  the  flow.  Thus, 
for  example,  on  the  X-24C-10D  the  side  and  bottom  of  the  vehicle  had  to  be 
handled  as  separate  panels  with  separate  surface  splines  to  obtain 
satisfactory  results. 

Streamlines  do  not  now  cross  panel  boundaries  because  data  must  be 
extrapolated  from  the  centroid  of  the  last  element  in  the  panel  to  the 
panel  edge,  then  as  the  first  centroid  in  the  new  panel  in  general  does 
not  fall  on  a  continuation  of  the  streamline,  the  data  must  be  interpolated 
in  the  circumferential  direction  and,  finally,  matched  with  the  data  saved 
on  the  next  set  of  centroids.  While  the  problem  is  not  insurmountable, 
the  risk  of  error  is  great. 

Finally,  a  flat  windward  surface  will  have  to  be  supplied  with  a  pressure 
distribution  from  some  other  source.  It  is  recommended  the  Mark  IV  HABP 
be  retained  to  estimate  heating  on  strakes  and  fins  where  other  codes  have 
proven  unsuccessful. 
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SECTION  III 

STUDIES  WITH  THE  DEJARNETTE  PROGRAM 
1.  PROGRAM  DESCRIPTION 

The  aerodynamic  heating  program  developed  by  Fred  R.  DeJarnette 
(References  4  and  9  through  1.4)  and  subsequently  modified  by  the  AFFDL 
(Reference  44)  is  described  in  this  section.  The  program  applies  the  small 
cross-flow  approximation  to  three-dimensional  boundary- layer  equations 
and  utilizes  the  axisymmetric  analog  technique  to  predict  laminar, 
transitional  and  turbulent  heating  rates  to  arbitrary,  blunt  nosed, 
three  dimensional  bodies  at  angle  of  attack  in  supersonic  and 
hypersonic  flow. 

Body  geometry  is  generated  by  a  bivariate,  cubic  spline  fit  of  the 
surface  coordinates.  For  simple  geometries  the  coordinates  are  calculated 
by  the  program.  For  complex  shapes  the  body  coordinates  can  be  read  into 
the  program.  Inviscid  streamlines  are  calculated  by  either  the  method 
of  steepest  descent  or  Euler's  momentum  equation.  The  surface  pressure 
is  either  calculated  by  the  modified  Newtonian  technique  or  read  into 
the  program  from  an  external  source.  Laminar  heating  rates  are  based  on 
the  locally  similar  boundary  layer  analysis  of  Cohen  (Reference  15)  and 
Beckwith  and  Cohen  (Reference  16) .  Turbulent  heating  rates  are  predicted 
by  applying  the  von  Kannan  form  of  the  Reynolds  analogy  to  the  semi- 
empirical  correlation  of  Spalding  and  Chi  (Reference  17)  using  the 
momentum  thickness  calculated  by  the  integral  method  of  Reshotko  and 
Tucker  (Reference  18) .  The  beginning  and  end  of  boundary  layer  transition 
is  specified  either  by  geometric  location,  momentum  thickness  Reynolds 
number  or  integrated,  local  unit  Reynolds  number  along  a  streamline. 
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Transitional  heating  rates  are  a  weighted  average  of  local  laminar 
and  turbulent  values  following  the  approach  of  Dhawan  and  Narasimha 
(Reference  19).  Entropy  layer  swallowing  effects  are  accounted  for  by 
a  modified  form  of  Maslen's  method  (Reference  20)  for  inviscid,  axisymmetric 
flows.  Either  perfect  or  equilibrium  air  properties  may  be  used. 

Better  understanding  of  the  program  can  be  achieved  if  its  various 
functions  are  described  separately. 

a.  Body  Description  -  Incorporated  in  the  program  is  a  bivariate, 
cubic  spline  fit  of  the  body  radius  as  a  function  of  the  longitudinal 
and  circumferential  coordinates  X  and  Body  coordinates  are  either 
calculated  by  or  read  into  the  program. 

A  special  provision  for  handling  body  data  originally  intended  for 
the  Supersonic-Hypersonic  Arbitrary  Body  Program  (Reference  1)  is  also  included. 

b.  Inviscid  Flow  Field  Calculations 

1.  Surface  Pressure  -  Surface  pressures  are  either  computed  by 
modified  Newtonian  theory  or  read  into  the  program.  For  80-degree  swept 
delta  wing  configurations,  a  pressure  correlation  formula  from  Reference.  21 
is  also  available.  Pressures  around  the  body  circumference  at  several 
axial  stations  are  used  to  generate  a  bivariate  cubic  spline  function  for 
the  pressure  distribution  at  any  point  on  the  body.  The  technique  is  similar 
to  that  used  to  describe  the  body  geometry.  In  "shadowed"  regions,  P  =  P® 

is  used. 

2.  Boundary-Layer-Edge  Gas  Properties  -  The  flow  field  calculation 
is  based  on  the  entropy  and  pressure  at  the  boundary  layer  edge.  The  entropy 
can  be  determined  by  either  of  two  different  techniques.  The  first  method 
simply  assumes  that  the  boundary-layer-edge  entropy  is  equal  to  the 
entropy  just  downstream  of  the  normal  shockwave.  This  assumption  is 
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invalidated  once  the  boundary  layer  has  swallowed  a  significant  portion 
of  the  entropy  layer.  The  second  method  incorporates  an  entropy  layer 
swallowing  technique  (described  in  detail  in  References  A  and  1A)  to 
determine  the  local  entropy.  The  usual  thin  boundary  layer  assumption 
sets  the  pressure  at  the  boundary  layer  edge  equal  to  the  surface  pressure. 

For  ideal  gases  surface  pressure  and  entropy  are  sufficient  to 
determine  all  the  other  gas  properties  using  isentropic  flow  relations. 

For  equilibrium  air  properties  the  correlation  formulas  of  Cohen 
(Reference  22)  are  used. 

3.  Streamline  and  Scale  Factor  Calculations  -  In  order  to  apply 
the  axisymmetric  analog  technique,  described  in  Section  1II-C,  to  three- 
dimensional  boundary-layer  equations,  it  is  necessary  to  compute,  the 
effective  body  radius  (scale  factor) .  This  requires  a  knowledge  of 
streamline  direction  and  curvature.  Two  differen  .  streamline  techniques 
are  available  in  the  program.  The  first  one  depends  on  the  body  orientation 
and  geometry  only  and  is  referred  to  as  the  method  of  steepest  descent 
in  Reference  23  (  and  throughout  the  remainder  of  this  report)  and  Newtonian 
streamlines  in  Reference  2A.  The  other  streamline  method  is  based  on  Euler's 
equation  written  in  streamline  coordinates.  This  latter  method  uses  both 
body  geometry  and  pressure  distribution  to  determine  the  streamline  geometry 
(Reference  11  and  13) . 

c.  Viscous  Boundary-Layer  Calculations  -  The  general,  three-dimensional 
boundary- layer  equations  are  reduced  to  equivalent  axisymmetric  equations 
by  using  the  small-croes-f low  assumption  of  Vaglio-Laurin  (Reference  25) . 

The  constraint  necessary  to  assure  that  the  three-dimensional  equations 
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can  be  reduced  to  axisymmetric  farm  in  that 


where  k  is  the  curvature  of  the  boundary  layer  edge  streamline  in  the 

plane  tangent  to  the  body  surface.  The  condition  |(u/ue)2  -  oe/pj  -*0 

l  j 

is  closely  approximated  when  the  wall  is  highly  cooled  and  ^/(Hq-H^P  <<  1. 
Experience  has  shown  that  this  last  constraint  may  be  overly  restrictive. 

The  axisymmetric  equations  that  are  obtained  can  be  solved  by  any 
method  applicable  to  a  body  of  revolution  at  zero  incidence.  The  distance 
along  the  streamline  is  interpreted  a3  the  distance  along  an  equivalent 
axisymmetric  body.  The  scale  factor  for  the  surface  coordinate  normal  to 
the  streamline  is  interpreted  as  the  radius  of  an  equivalent  axisymmetric 
body  at  zero  incidence.  This  method  is  commonly  referred  to  as  the 
axisymmetric  analog  technique. 

d.  Heat  Transfer  Rate  Predictions 

1.  Stagnation  Region  Heat  Transfer  Rate  -  Using  the  small  cross 
flow  approximation,  the  heating  rate  along  an  inviscid  surface  streamline 
is  obtained  from  Lees'  heat  transfer  equation  for  a  cold  wall  (Reference  26) 
on  an  equivalent,  axisymmetric  body  at  zero  incidence  by  replacing  the 
body  radius  by  the  scale  factor  and  the  distance  along  the  body  surface 
by  the  distance  along  the  streamline.  The  limit  of  this  equation  is  taken 
as  the  stagnation  point  is  approached.  Finally,  the  resulting  equation 
is  modified  so  that  it  produces  results  compatible  with  experimental 
results  in  the  limit  of  two-dimensional  and  axisymmetric  stagnation 
points  (Reference  11). 
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2.  Surface  Heating  Rates 

(1)  Laminar  Boundary  Layers 

Laminar  heating  rate  ratios  are  determined  by  applying 
the  axisymmetric  analog  to  solutions  of  the  locally  similar  boundary 
layer  equations  for  relatively  cold  walls  given  by  Beckwith  and  Cohen 
in  Reference  16. 

(2)  Turbulent  Boundary  Layers 

The  axisymmetric  analog  is  applied  to  a  modified  form 
of  the  integral  method  of  Reshotko  and  Tucker  (Reference  18)  to  obtain 
the  turbulent  momentum  thickness.  Then  the  momentum  thickness  is  used 
to  calculate  the  local  skin  friction  coefficient  from  the  semi-empirical 
correlation  of  Spalding  and  Chi  (Reference  17).  Finally,  skin  friction 
coefficients  are  converted  to  turbulent  heating  rates  through  the 
von  Karman  form  of  Reynolds  analogy. 

(3)  Transition  Region  Heating  Rates 

The  transition  region  heating  rates  are  based  on  a 
weighted  average  of  the  local  laminar  and  turbulent  heating  rates.  The 
method  (based  on  the  results  of  Dhawan  and  Narasimha  in  Reference  19) 
uses  a  Gaussian  distribution  for  the  weighting  function  which  is  based 
on  the  geometric  locations  of  the  beginning  and  end  of  transition.  The 
analysis  does  not  predict  transition  ensjet  or  extent.  Instead,  the 
beginning  and  end  of  transition  is  specified  in  the  input  data  by  one  of 
three  options  which  are:  geometric  location,  value  of  the  local  momentum 
thickness  Reynolds  number,  or  value  of  the  integrated  local  unit  Reynolds 


number 
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2.  PROGRAM  LIMITATIONS  AND  SOURCES  OF  ERROR 

Since  the  body  geometry  is  approximated  by  a  bivariate  cubic  spline 
fit,  there  exist  certain  restrictions  that  pertain  to  the  input  geometry. 

One  restriction  is  that  the  body  radius  must  be  single  valued.  This  means 
that  a  cross  section  such  as  shown  in  the  sketch  below  cannot  be  computed 
because  the  body  radius,  f,  is  multi-valued  on  certain  portions  of  the  body. 


In  addition,  the  program  is  restricted  to  relatively  smooth  aerodynamic 
configurations  such  as  shown  in  Figure  5.  In  general,  only  blunt  nose 
configurations  can  be  computed  although  there  is  an  input  option  that 
sets  the  inviscid  streamline  pressure  equal  to  the  input  surface  pressure. 
With  this  option  and  the  Euler  streamline  option,  reasonably  good  results 
can  be  obtained  for  sharp  nose  bodies. 

In  theory,  the  bivariate  cubic  spline  function  can  be  used  to  describe 
very  general  body  shapes.  However,  when  a  cylindrical  coordinate  system  is 
used,  as  in  the  present  program,  computations  are  inaccurate  on  planar 
type  bodies.  This  is  due  to  the  non-linearity  of  the  second  derivative 
of  the  square  of  the  body  radius  with  respect  to  the  circumferential  angle. 
The  result  is  that,  on  the  flat  section,  oscillations  in  the  body  surface 
slope  occur  that  induce  errors  in  the  streamline  solution. 

When  reading  body  geometry  data  into  the  program,  the  axial  stations 
are  not  required  to  be  equally  spaced.  However,  the  spline  fit  in  the  axial 
direction  is  generally  more  accurate  when  the  axial  stations  are  spaced 
approximately  equal.  Although  the  circumferential  body  points  do  not  have 
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to  be  equally  spaced,  it  is  required  that  the  first  body  point  be  on  the 
windward  streamline  (<J>  =  0°)  followed  by  other  points  located  at 
monotonically  increasing  values  of  $  until  the  final  body  point  is 
reached  on  the  leeward  streamline  (4>  =  180°). 

Computer  storage  restrictions  limit  the  number  of  body  points  that 
can  be  input  to  the  program.  Only  a  maximum  of  20  cross  sections,  each 
having  a  maximum  of  20  circumferential  body  points,  is  allowed.  In  addition, 
no  restart  capability  exists.  Therefore,  the  size  and  shape  of  aerodynamic 
configurations  that  can  be  considered  is  somewhat  restricted. 

The  inviscid  flow  field  properties  at  the  edge  of  the  boundary  layer 
are  computed  from  the  local  entropy  and  pressure.  Therefore,  the  flow- 
field  accuracy  is  dependent  on  the  accuracy  of  the  surface  pressure 
distribution  and  the  validity  of  the  entropy  swallowing  technique. 

Pressure  data  can  either  be  read  into  the  program  or  be  calculated  from 
modified  Newtonian  theory.  In  general,  pressures  on  blunt  slender  bodies 
are  not  adequately  described  by  the  Newtonian  assumption  nor  are  pressures 
over  locally  flat  portions  of  a  three-dimensional  body.  In  addition,  leeside 
pressures  are  usually  quite  different  from  those  obtained  from  a  Newtonian 
calculation.  If  accurate  pressure  distributions  are  available,  the 
calculated  inviscid  flow  field  characteristics  and  streamline  orientation 
should  closely  approximate  the  physical  situation. 

Another  source  of  error  can  develop  whenever  a  non-monotonic  pressure 
distribution  (in  either  the  X  or  <j>  coordinate  directions)  is  used.  It  has 
been  observed  that  whenever  a  non-mono tonic  pressure  distribution  is  used. 
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the  scale  factor  behaves  strangely  downstream  of  the  pressure  distribution 
inflection  point. 

Streamlines  computed  from  the  method  of  steepest  descent  are  based 
solely  on  a  geometric  concept  and  cannot  be  expected  to  be  physically 
representative  when  overriding  aerodynamic  forces  are  present.  Since 
the  steepest  descent  streamline  calculation  is  independent  of  surface 
pressures,  its  accuracy  for  axisymmetric  bodies  at  zero  degrees  angle 
of  attack  cannot  be  degraded  by  an  inaccurate  pressure  distribution. 
However,  heating  rates  along  these  streamlines  on  axisymmetric  bodies 
at  zero  degrees  angle  of  attack  can  be  in  error  due  to  inaccurate 
inviscid  flow  conditions  calculated  from  an  inaccurate  pressure 
distribution. 

Streamlines  determined  from  Euler’s  equation  are  dependent  on  the 
pressure  distribution  as  well  as  the  body  geometry.  Consequently,  any 
error  in  the  presoure  field  will  affect  streamline  direction  and 
curvature,  because  the  equivalent  body  radius  (scale  factor)  is  dependent 
solely  on  the  position  and  curvature  of  the  streamline,  an  accurate 
streamline  pattern  is  essential  for  the  success  of  the  program. 

The  accuracy  of  the  axisymmetric  analog  method  depends  on  the  degree 
to  which  the  governing  assumptions  of  small  cross  flow  and/or  cool  wall 
temperatures  are  satisfied. 

While  the  program  accounts  for  the  effect  of  an  adverse  surface- 
pressure  gradient,  no  provision  is  made  to  account  for  the  effect  of  a 
wall  temperature  gradient  for  either  laminar  or  turbulent  boundary  layers. 
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The  turbulent  boundary  layer  analysis  is  expressed  in  terms  of  the 
boundary  layer  form  factor  and  momentum  thickness.  Development  of  the 
analysis  along  these  lines  is  due  solely  to  the  need  for  obtaining  a 
reasonably  compact  engineering  technique  for  solving  the  boundary 
layer  equations. 

3.  SAMPLE  TEST  CASES 

In  this  section  experimental  heat  transfer  data  are  compared  with 
theoretical  predictions  for  four  different  test  models.  Program  accuracy 
and  some  of  its  limitations  are  noted  where  appropriate.  A  summary  of 
each  body  geometry,  its  attitude  and  the  corresponding  flow  conditions 
is  given  in  Table  I. 

1.  25°  Half-Angle  Cone  at  a  -  0° 

This  is  a  very  simple  case  intended  to  demonstrate  the  accuracy 
of  the  program  for  the  aero  cross-flow  situation.  It  has  a  spherically 
blunted  nose  with  a  1/ 2-inch  nose  radius.  Freestream  Mach  number  was 
7.95  and  the  freestream  unit  Reynolds  number  was  3.96  x  10^  per  foot. 
Detailed  freestream  flow  conditions  are  given  in  Table  I.  The 
experimental  heat  transfer  coefficients  are  based  on  a  recovery  factor 
of  0.85,  where  the  recovery  factor  is  defined  as  (Taw  -  Te)/(To  -  Te) . 

For  the  calculation  of  local  flow  conditions,  an  isentropic  expansion 
From  the  model  stagnation  point  to  the  local  experimental  pressure  level 
was  assumed.  Test  data  were,  reported  by  Bushnell  et  al,  in  Reference  27, 
Both  analytical  and  experimental  heat-transfer  coefficients  were  ratioed 
to  the  stagnation  point  value  given  in  Reference  27  which  is  h,^  “ 

34  x  10“3  (Btu/f t^-sec-°R) . 

In  Figure  6,  a  comparison  of  experimental  surface  pressures  ratioed 
to  stagnation  point  pressure  is  made  with  a  spline  fit  of  the  modified 
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Newtonian  pressure  distribution.  Clearly,  the  modified  Newtonian  result 
is  in  considerable  error  as  is  to  be  expected.  Also  shown  in  Figure  6 
are  pressures  determined  analytically  by  the  method  described  in 
Reference  28.  The  spline  fit  of  both  pressure  distributions  produces 
a  pressure  dip  near  S/R  =  1.5.  This  is  due  to  limitations  inherent  in 
the  bivariate  cubic  spline  fitting  technique. 

Using  the  fitted  Newtonian  pressure  distribution,  the  calculated 
heating  rates  are  essentially  identical  for  streamlines  calculated 
either  by  the  steepest  descent  method  or  by  the  Euler  equation  method. 
Results  of  both  methods  are  shown  as  the  broken  line  in  Figure  7. 

Because  the  Newtonian  pressure  approximation  does  not  account  for  the  over 
expansion  region  just  downstream  of  the  sphere-cone  junction,  the 
predicted  heating  rates  between  S/R  *  1.2  and  S/R  =  7.0  are  noticeably 
higher  than  the  data.  Using  the  spline  fit  of  the  pressures  from 
Reference  28,  somewhat  different  heating  rates,  shown  in  Figure  7  as  the 
dashed  line,  result.  The  accuracy  of  the  predicted  heat  transfer 
distribution  is  contingent  on  the  accuracy  of  the  surface  pressure 
distribution.  Consequently,  use  of  the  more  accurate  pressure  field 
results  in  improved  agreement  between  the  experimental  and  calculated 
heating  rate  distribution  as  can  be  seen  in  Figure  7.  However,  results 
given  by  the  present  program  differ  somewhat  from  the  results  of  Bushnell, 
et  al.  In  Reference  27.  This  is  due  to  the  inaccuracy  of  the  spline  fit 
of  the  pressure  distribution  between  S/R  =  1.0  and  S/R  =  2.8,  shown  in 
Figure  6  as  a  dotted  line. 
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2.  0. 7:1.0  Elliptic  Cone 

The  second  test  case  considered  is  that  of  a  blunt,  elliptical 
cone  having  a  semi -minor  to  semi-major  axis  ratio  of  0,7.  Data  for 
this  case  were  taken  from  experimental  results  reported  by  Hillsamer 
and  Rhudy  in  Reference  29.  A  description  of  the  geometric  properties 
of  the  test  vehicle  is  given  in  Figure  8.  Nominal  freestream  flow 
conditions  are  given  in  Table  I.  Experimental  and  analytical  heat- 
transfer  coefficients,  h,  are  based  on  total  temperature, 

h  -  qw/'(T0  -  Tw).  (2) 

Nominal  values  are  To  “  1910  (°R)  and  Tw  ■  470  (°R) .  Local  flow  field 
calculations  utilize  the  entropy  layer  swallowing  process  in  conjunction 
with  a  modified  Newtonian  pressure  distribution. 

Heating  rates  on  the  windward  streamline  <4>  ■  0°)  are  shown  in 
Figure  9  for  a  »  0°  and  oi  -  30°  angle  of  attack.  At  a  °  0°  it  can  be 
seen  that  neither  the  steepest  descent  nor  the  Euler  streamline 
computations  give  very  accurate  heating  rates  when  compared  with  experimental 
data.  As  in  the  25°  cone  case,  the  inaccuracy  is  due  primarily  to 
limitations  of  the  modified  Newtonian  pressure  distribution  at  small 
inclination  angles.  However,  the  steepest  descent  method  provides  more 
accurate  heating  rates  than  the  Euler  equation  method.  This  is  because 
the  effect  of  an  inaccurate  pressure  field  1b  compounded  when  using  the 
Euler  equation  to  determine  streamline  orientation  and  body  scale  factor. 

In  the  heat-transfer  calculation  not  only  are  the  inviscid  flow  values 
in  error,  but  so  is  the  integration  of  the  Euler  equation  which  is  dependent 
on  the  inaccurate  pressure  distribution.  By  comparison,  inaccuracy  in 
the  steepest  descent  method  (for  this  case  at  $  “  0°)  is  due  entirely  to 
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erroneous  flow  field  values  calculated  from  the  spline  fit  of  the 
modified  Newtonian  pressure  distribution. 

The  circumferential  calculations  at  a  *  0°  are  shown  in  Figure  10. 

It  can  be  seen  that  the  heating  distribution  approaches  the  experimental 
values  near  <J)  «  90°  where  the  body  surface  is  more  steeply  inclined  to 
the  freestream  flow  direction.  Agreement  between  data  and  the  steepest 
descent  theory  is  seen  to  improve  with  increasing  distance  downstream  of 
the  stagnation  point.  At  station  14.5  the  agreement  is  seen  to  be  quite 
good  especially  near  <f>  -  90°.  As  for  heating  distributions  calculated 
from  the  Euler  equation,  the  agreement  with  data  is  uniformly  poor. 

At  a  ■  30°  angle  of  attack,  both  heating  rate  prediction  techniques 
give  values  that  agree  quite  well  with  experiment.  In  Figure  9  the 
agreement  between  data  and  the  method  of  steepest  descent  is  seen  to  be 
slightly  better  than  the  method  based  on  Euler's  equation.  The  same 
results  can  also  be  seen  in  Figure  10.  The  agreement  between  data  and 
theory  on  the  leeside  of  the  model  is  seen  to  be  relatively  poor  as  would 
be  expected. 

3.  Blunt,  80°  Delta  Wing 

This  test  case  is  a  blunt,  80°  delta  wing  at  a  »  40°  angle  of  attack, 
Mw  -  7.97  and  unit  freestream  Reynolds  number  Re,,,  ■  1.94  x  10^  per  foot 
(Reference  21).  Detailed  freestream  flow  conditions  are  given  in  Table  I. 

The  model  is  29.4  inches  long  and  has  a  1/8  inch  nose  radius.  This  test 
case  demonstrates  the  deleterious  effect  that  a  flat  surface  has  on  heat- 
transfer  calculations. 

The  experimental  heat-transfer  coefficients  shown  in  Figures  11  and  12 
were  determined  from  the  semi-infinite  fl/'t  plate  heat  conduction  equation. 
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Since  all  heat  transfer  data  were  taken  at  X/R  >  40,  it  was  assumed 
that  the  blunt  nose  induced  entropy  layer  had  been  completely  swallowed 
by  the  boundary  layer  upstream  of  the  S/R  =  40  location.  Consequently, 
experimental  heat-transfer  coefficients  were  based  on  local  flow 
conditions  calculated  from  experimental  surface  pressures  in  conjunction 
with  oblique  shockwave  theory. 

Numerical  heat-transfer  calculations  were  made  using  the  entropy 
layer  swallowing  analysis  and  the  experimentally  determined  pressure 
distribution.  All  calculations  used  a  perfect  gas  with  Y  *  1.4. 

In  Figure  11  it  can  be  seen  that  on  the  windward  streamline  at 

a  “  40°  there  is  fairly  good  agreement  between  analysis  and  experiment. 

However,  in  Figure  12  the  agreement  between  theory  and  data  at  spanwise 

locations  is  quite  poor  in  some  instances.  The  worst  agreement  occurs 

between  Z/Z  “  0.12  and  Z/Z  =  0.4  at  X  “12.  In  this  region  the 
max  max  vn 

body  radius  and  its  first  and  second  derivatives  with  respect  to  X  and 
<f>  were  consistently  more  in  error  than  in  nearby  regions  where  the 
agreement  between  fitted  and  sxact  results  was  better.  Therefore,  it 
is  suspected  that  the  poor  agreement  between  data  and  the  calculated 
results  is  due  to  the  inaccurate  surface  fit.  The  primary  drawback 
regarding  cubic  splines  is  that  the  second  derivative  of  the  function 
being  fitted  must  vary  linearly  with  respect  to  the  independent  variable. 
When  describing  a  flat  section  using  a  cylindrical  coordinate  system, 
the  second  derivative  of  the  square  of  the  body  radius  varies  non-linearly 
with  the  circumferential  angle.  ThuB,  the  program  is  limited  to 
non-planar  type  bodies.  On  the  flat  section,  oscillations  in  the  body 
surface  slope  occur  that  induce  errors  in  the  streamline  solution  for 
the  scale  factor. 
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4.  10°  Half-Angle  Cone 

The  last  test  case  is  a  blunt,  10°  half-angle  cone  at  a  =  50° 
angle  of  attack,  =  8.0  and  unit  freestream  Reynolds  number,  Reoa  of 
3.7  x  10^  per  foot  (Reference  30).  Detailed  freestream  flow  conditions 
are  given  in  Table  I.  The  nose  is  a  1.0-inch-radius  sphere,  and  the  model 
is  25.24  inches  long. 

This  test  case  was  chosen  to  demonstrate  the  transition  region  heating 
rate  prediction  capabilities  of  the  program  for  a  simple  body  shape. 

The  experimental  heat-transfer  coefficients  shown  in  Figures  13  and  14 
(Reference  30)  were  determined  from  the  semi- inf ini te  flat-plate  heat 
conduction  equation  with  the  simplifying  assumption  that  Taw  =  TQ. 
Theoretical  results  were  adjusted  for  an  adiabatic  wall  temperature  equal 
to  the  total  temperature.  Theoretical  stagnation-point  heating  rates 
were  set  equal  to  experimental  results  so  that  a  comparison  with  data  could 
be  made.  A  perfect  gas  withy  “1.4  was  used  for  both  the  analytical  and 
experimental  calculations. 

In  the  program  transitional  heat-transfer  rates  are  calculated  as  a 
weighted  average  of  the  local  laminar  and  turbulent  heating  rates.  The 
transition  region  heating  rate  is  defined  as 

*w-Vu„  +  "f  «wturt'  w-  <3) 

Calculations  for  the  weighting  factor,  Wf,  are  based  on  the  flat-plate, 
low-speed,  semi-empirical  analysis  of  Dhawan  and  Narasimha  (Reference  19). 
This  analysis  uses  a  Gaussian  distribution  for  Wf,  which  is  given  by 
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DeJarnette  (Reference  11)  assumes  that  this  distribution  can  be  applied 

directly  to  compressible  flows  having  streamline  divergence.  This  is 

tantamount  to  assuming  that  the  inviscid  flow  conditions  at  the  boundary 

layer  edge  do  not  vary  from  the  beginning  to  end  of  transition,  and  that 

the  streamline  divergence  is  constant  over  the  transition  region. 

DeJarnette  further  assumes  that  the  distribution  for  can  also  be 

accurately  determined  if  either  local  momentum  thickness  Reynolds  number, 

s  Pe  ue 

Pe  ue  Qm^he*  or  local  integrated  unit  Reynolds  number,  /  -  ds, 

o 

are  substituted  for  the  physical  X-locations  in  Equation  (4).  For  such 
a  substitution  to  be  valid,  both  momentum  thickness  and  integrated  unit 
Reynolds  number  must  vary  linearly  with  X. 

All  three  program  options  for  predicting  transition  region  heating 
rates  were  used  for  this  case.  Results  using  the  geometric  location 
option  and  the  integrated  local  unit  Reynolds  number  option  were  essentially 
identical  and  are  shown  in  Figure  13.  Reasonably  good  agreement  between 
theory  and  experiment  can  be  seen.  Results  from  the  momentum  thickness 
Reynolds  number  option,  also  shown  in  Figure  13,  differed  considerably 
from  data  and  results  from  the  other  two  options.  The  reason  for  this 
is  that  program  calculations  for  the  local  momentum  thickness  Reynolds 
number  showed  a  very  non-linear  variation  with  X.  By  comparison,  the 
local  integrated  unit  Reynolds  number  varied  (in  this  case  at  least) 
linearly  with  X.  Examination  of  the  momentum  integral  equation  for  high¬ 
speed  flow  over  axisymmetric  bodies,  Equation  (5)  below. 
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reveals  that  momentum  thickness,  0m!,  is  a  complicated  function  of  flow 
conditions  and  body  geometry.  Even  for  the  simple  case  of  laminar 
flow  over  a  flat  plate,  where 

dQm 

a  l/^/xT  (6) 

the  momentum  thickness  does  not  vary  linearly  with  X.  Within  the 
transition  region  the  relation  between  momentum  thickness  and  X  is 
certainly  not  linear.  We  can  conclude  then,  that  the  momentum  thickness 
Reynolds  number  option  is  a  poor  choice  to  use  for  predicting  heating 
rates  in  the  transition  region. 

Calculated  apanwise  heating  rates,  shown  in  Figure  14,  agree  quite 
well  with  experimental  data  taken  in  the  transition  region.  The  local 
integrated  unit  Reynolds  number  option  was  used  to  define  the  beginning 
and  end  of  transition.  End-point  values  were  based  on  windward  streamline 
data.  At  apanwise  locations  laminar  theory  compares  less  favorably 
with  data  than  do  transition  region  results.  In  turbulent  flow  the 
effects  of  streamline  divergence  are  much  less  than  for  laminar  flow. 

This  effect  is  carried  over  into  the  transition  region  to  give  good 
agreement  between  analytical  and  experimental  results  at  X/L  »  0.7 
and  X/L  »  0.9  as  shown  in  Figure  14. 

One  of  the  inherent  drawbacks  of  a  boundary  layer  integration  along 
a  streamline  is  that  streamlines  tend  to  spread  or  converge  depending  on 
body  shape  and  pressure  distribution.  There  is  no  way  of  knowing 
beforehand  just  which  stagnation  region  streamlines  should  be  chosen  in 
order  to  provide  a  uniform  circumferential  streamline  distribution  at  a 
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certain  body  cross  section  downstream  of  the  stagnation  region.  Also, 
any  set  of  stagnation  region  streamlines  that  produce  a  uniform 
circumferential  distribution  of  streamlines  at  one  body  station  will 
probably  not  produce  a  very  good  circumferential  streamline  distribution 
at  a  significantly  different  body  station  if  sizeable  spanwise  surface 
or  pressure  gradients  are  present. 
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SECTION  IV 

THREE-DIMENSIONAL  SHOCK-CAPTURING  TECHNIQUE 
1.  BACKGROUND 

While  empirical  techniques  have  been  very  successful  in  predicting 
the  surface  pressure  on  large  areas,  the  increased  emphasis  on  heating 
rates  requires  better  streamline  definition  which,  in  turn,  requires 
better  pressure  detail.  Some  more  nearly  "exact"  program  is  required 
to  supply  such  pressure  detail.  It  should  be  noted  that  thi3  is  not 
quite  as  difficult  as  might  appear,  because  surface  pressure  is  relatively 
insensitive  to  errors  in  the  flow-field  calculation  (one  of  the 
reasons  for  the  success  of  empirical  approximations).  The  Northrop 
Three-Dimensional  Method  of  Characteristics  (Reference  31)  has  been 
used  for  this  purpose,  although  there  are  geometry  limits  and  the 
development  of  a  strong  cross-flow  shock  will  result  in  failure  of  the 
characteristic  code.  A  three-dimensional  time-dependent  solution  is  much 
too  costly  and  the  Grumman  STEIN  code  (a  three-dimensional  shock- tracking 
technique)  (Reference  32)  was  not  locally  available  and  was  reported  to 
be  very  large  and  slow  (Reference  33).  The  Three  Dimensional  Shock 
Capturing  Technique  (SCT)  developed  by  Solomon  (Reference  6)  is 
available.  The  SCT  calculates  the  flow  in  a  supersonic  region  by  solving 
the  steady,  inviscid  equations  written  in  a  weak  conservation  form.  With 
the  equations  in  this  form,  imbedded  shocks  are  located  in  an  approximate 
manner;  the  bow  shock  is  explicitly  located  and  tracked  by  the  program. 

The  SCT  code  is  able  to  handle  steep  entropy  gradients  near  the  body, 
vehicles  which  are  not  symmetric  about  the  pitch  plane  and  perform  real- 
gas  as  well  as  perfect-gas  calculations. 
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Because  the  program  employs  an  explicit  finite  differencing  scheme 
which  advances  the  solution  along  the  longitudinal  axis,  the  storage 
requirements  are  minimized.  As  each  surface  is  completed,  it  becomes 
the  initial  value  surface  for  the  next  step,  thus  the  program  may  be 
halted  at  any  station,  the  surface  saved  and  the  program  restarted  from 
that  location. 

All  geometry  calculations  are  carried  out  in  a  subroutine  "BODY" 
which  computes  body  radius  as  a  function  of  meridional  angle;  the  first 
and  second  derivatives  of  body  radius  with  respect  to  meridional  angle; 
the  first  and  second  derivatives  with  respect  to  the  longitudinal  direction 
and  the  cross  derivative  with  respect  to  longitudinal  direction  and  meridional 
angle.  All  other  required  manipulation  of  geometry  data  is  in  a  separate, 
permanent  subroutine.  Thus,  if  a  body  is  required  which  cannot  be 
described  by  the  existing  routine,  a  new  BODY  subroutine  can  be  used  in 
place  of  the  original  routine. 

As  delivered,  the  BODY  routine  was  designed  to  calculate  a  multiconic 
body  with  cuts,  flaps,  and  flares.  This  routine  was  used  for  the  blunt 
and  sharp  biconic  studies  described  in  the  next  section. 

For  the  type  aerodynamic  study  of  interest  here,  however,  different 
bodies  were  required.  The  first  non-conical  shape  attempted  was  a  blunted 
ogive-cylinder  for  which  extensive  experimental  data  was  available 
(References  34  and  36) .  The  new  subroutine  is  described  in  Appendix  A. 

Some  results  of  these  calculations  are  discussed  later  in  this  section. 

More  complex  lifting  configurations  required  a  more  elaborate  program. 
Details  of  this  code  are  included  in  Appendix  B  of  this  report.  This 
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subroutine  calculates  the  body  geometry  for  any  body  which  can  be  described 
by  a  set  of  circular  arcs  and  straight  lines.  Elliptic  arcs  could  also  be 
included  but  have  not  been  as  yet.  The  subroutine  assumes  the  slope  is 
continuous  at  the  junction  of  each  segment  in  the  meridional  direction; 
the  body  may  be  discontinuous  in  the  longitudinal  direction. 

The  SCT  is  a  fast,  efficient  code  and  the  comparisons  with  experimental 
data  are  excellent.  The  code  is  not  completely  automatic;  some  knowledge 
of  the  flow  behavior  and  behavior  of  the  equations  is  necessary  to  tailor 
the  run  set  up,  but  it  has  failed  only  when  a  compression  forced  the 
axial  flow  to  go  subsonic  or  downstream  of  large  expansions  in  either  the 
longitudinal  or  circumferential  directions  which  caused  numerical 
instabilities  in  the  solution.  Angles  of  attack  in  excess  of  30°  for  15° 
half-angle  cones  have  been  calculated  and  a  few  demonstration  cases  with 
a  small  yaw  angle  were  completed. 

2.  CONE  STUDY 

The  SCT  Program  was  developed  primarily  tc  study  three-dimensional 
flow  over  conical  bodies  with  flats,  flaps  and  flares.  Therefore,  a 
series  of  Liconic  shapes  were  selected  as  the  initial  cases  to  gain  some 
experience  with  the  program  before  attempting  to  extend  the  program's 
capabilities.  This  study  used  a  5°/ 10°  biconic  with  nose  bluntness 
ratios  (Rn/RBase)  from  2 %  to  20%  and  angles  of  attack  up  to  25°.  The 
experimental  study  was  on  stability  and  no  heat  transfer  comparisons 
were  available.  Results  of  this  study  were  presented  at  the  AIAA  4th 
Atmospheric  Flight  Mechanics  Conference  (Reference  35) . 

One  computational  problem  is  worth  discussing.  At  small  angles  of 
attack  (less  than  3°),  it  was  possible  to  complete  the  full  solution  on 
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the  computer  without  intervention;  for  larger  angles  it  was  necessary  to 
stop  calculation  about  three  nose  radii  beyond  the  biconic  junction  and 
restart  the  solution  at  that  point.  The  reason  for  this  is  the  manner 
in  which  the  program  computes  wall  boundary  conditions. 

It  is  possible  to  derive  two  equivalent  equations  for  the  derivative 
of  pressure  in  the  axial  direction  depending  on  which  substitutions  are 
made  in  the  basic  equation.  These  two  equations,  however,  may  result  in 
significant  differences  in  the  computed  pressures.  More  significantly, 
one  may  give  a  solution  in  regions  where  the  other  will  not.  The  two 
solutions  are  referred  to  as  "MOD  3"  and  "MOD  0"  solutions;  the  MOD  3 
solution  has  proven  to  be  the  most  dependable  for  bodies  of  revolution 
except  in  regions  of  surface  discontinuity  in  the  axial  direction. 

In  this  case,  the  MOD  3  solution  often  fails;  the  MOD  0  solution  gives 
good  solutions  in  this  case.  However,  the  MOD  0  solution  does  not,  in 
general,  compute  the  leeside  of  cones  at  large  angles  of  attack  successfully. 
Because  of  the  behavior  of  these  two  solutions,  the  default  value  of  the 
program  is  MOD  3  until  an  axial  discontinuity  is  reached,  i.e.»  Subroutine 
JUMP  is  called,  at  which  time  MOD  0  is  selected.  There  is  no  logic  to 
permit  the  user  to  control  the  selection  of  these  two  solutions. 

In  his  report,  Solomon  (Reference  6)  makes  no  effort  to  explain  why 
one  boundary  condition  works  when  the  other  does  not;  the  complexity  of 
the  equation  is  too  great  to  make  any  definitive  statements  here.  However, 
some  clarifying  generalities  can  be  made.  In  the  primary  pressure  method, 
MOD  3,  numerical  differences  are  required  for 
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p,  v£_,  uv,  vw,  and  v 
RB  Rb  Rb  rB 

while  in  the  alternate  boundary  condition,  MOD  0,  numerical  differences 

are  required  only  for 

P,  v  . 
w 

In  regions  of  sudden  large  expansion  MOD  3  breaks  down,  probably 
because  the  products  v^,  uv,  vw  display  very  steep  gradients  in  such  a 
region.  Dividing  the  circumferential  velocity  by  the  longitudinal  velocity 
results  in  a  parameter  which  is  not  changing  rapidly,  permitting  successful 
differences  to  be  taken  in  the  tangential  direction.  On  the  leeside 
of  cones  at  large  incidence  the  ratio  of  v/w  tends  toward  unity  because  of 
the  very  high  cross  flow  velocity  resulting  in  poor  differences  in  the 
tangential  direction. 

Using  the  above  guidelines,  it  would  seem  wise  to  force  a  MOD  0  type 
solution  any  time  the  ratio  v/w  is  small  (cross  flow  is  small  compared  to 
the  axial  flow  velocity) ,  and  p,  u  and  v  show  high  gradients  and 
force  MOD  3  solutions  any  time  the  cross  flow  approaches  the  axial 
velocity  (v/w  tends  toward  unity) . 

3.  OGIVE-CYLINDER  STUDY 

AFFDL  has  conducted  an  extensive  series  of  tests  at  AEDC  on  a  modular 
missile  configuration  (Reference  34)  to  study  fin  interactions.  Phase  III 
of  the  series  (Reference  36)  was  planned  to  supply  data  for  comparison 
with  the  SCT  Program.  A  surface  oil  flow  study  was  conducted  in  an  effort 
to  determine  streamline  locations  along  the  body.  Oil  flow  data  was  taken 
at  angles  of  attack  of  4°,  6°,  8°  and  10°.  The  fins  shown  in  Figure  15 
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were  removed  and  a  pitot  pressure  rake  and  a  total  temperature  probe 
were  installed  180°  apart  and  the  run  sequence  was  repeated.  Static 
pressures  on  the  surface  of  the  model  at  body  station  41  were  also  taken 
during  this  phase  of  the  test. 

A  series  of  ogive-cylinder  computer  runs  using  the  NSWC  code  with 
the  subroutine  described  in  Appendix  A  was  made  at  angles  of  attack  of  0°, 

2°,  4°,  5°,  6°,  7°,  8°,  9°,  10°  and  12°.  Computer  runs  using  the 
DeJarnette  heat  transfer  code  were  made  for  comparison  with  heat  transfer 
measurements  taken  in  other  phases  of  the  study.  No  effort  was  made  to 
calculate  flow  in  the  fin  region. 

Surface  static  pressures,  normalized  by  frees tream  static  pressure, 
are  shown  in  Figure  16.  The  SCT  pressures  compare  well  with  the  experimental 
values  on  the  windward  portion  of  the  model;  the  leeside  pressures  are 
not  well  matched,  nor  should  they  have  been  expected  to  match  due  to  the 
reasons  cited  in  Section  IV-2. 

Figure  17  shows  the  results  of  the  pitot  rake  survey  of  the  inviscid 
shock  layer.  The  program  predicts  the  pitot  pressure  successfully  in  the 
(J^O  and  90°  planes;  however,  viscous  effects  at  <j>  “  180°  are  responsible 
for  the  large  discrepancies  between  experimental  and  calculated  pressures 
close  to  the  body.  There  are,  however,  very  favorable  comparisons  near 
the  shock,  even  on  the  leeside.  The  inviscid  stream  angles  predicted  by 
the  SCT  were  consistently  lower  than  the  stream  angles  indicated  by  the 
oil  flow  experiment  as  expected.  Consequently,  the  inviscid  stream  angles 
were  used  in  the  correlation  theory  for  which  this  study  was  conducted 
and  showed  better  correlation  than  the  oil  flow  angles.  Additional  details 
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of  the  study  were  presented  at  the  AIM  15th  Aerospace  Science  Meeting 
(Reference  37) . 

4.  LIFTING  BODY  STUDY 

The  original  body  geometry  portion  of  the  SCT  permitted  only 
axisymmetric  bodies  with  cuts  and  flaps.  The  extension  to  ogive- 
cylinder  bodies  required  a  new  subroutine  to  describe  the  geometry,  but 
the  shape  was  still  axisymmetric.  However,  a  large  class  of  vehicles, 
such  as  lifting  bodies,  have  only  one  plane  of  symmetry  and  are  distinguished 
by  large,  flat,  or  nearly  flat  surfaces.  To  utilize  the  SCT  for  these 
bodies  a  new  geometry  subroutine  was  required. 

This  new  geometry  subroutine  must  permit  the  input  of  shapes  which 
vary  in  both  the  circumferential  and  longitudinal  directions.  The  data 
required  by  the  new  subroutine  should  be  as  compatible  with  the 
information  supplied  on  engineering  drawings  as  possible.  The  subroutine 
described  in  Appendix  B  was  developed  with  reference  to  the  X--24C-10  model 
drawings  which  supply  lofting  type  information. 

There  are  at  least  two  approaches  to  producing  the  required  body 
information.  The  obvious  one  is  an  analytic  description  of  the  geometry 
in  an  appropriate  coordinate  system,  the  other  involves  surface  fitting 
the  geometry  with  polynomials.  The  latter  approach  is  the  one  taken  by 
Marconi  (Reference  32).  In  such  a  scheme,  the  physical  body  must  be  mapped 
into  a  "mapped  space"  to  define  the  body,  then  into  the  computational  space. 

Although  there  is  a  great  deal  of  merit  to  the  surface  fitting  approach, 
the  purely  analytic  one  was  chosen  for  this  program.  The  basic  reason  for 
this  selection  was  the  ease  of  taking  loft-line  information  from  engineering 
drawings  and  using  them  in  a  computer  code.  In  the  pure  sense  of  the  term, 

33 


AFFDL-TR- 7 9-3001 


this  code  also  maps  from  cartesian  coordinates  to  cylindrical  coordinates 
to  a  rectangular  computational  coordinate  system;  the  most  significant 
difference  is  that  in  this  scheme  a  1:1  mapping  from  cartesian  to  cylindrical 
coordinates  is  always  used.  If  surface  fitting  is  the  objective,  the 
mapping  is,  in  general,  not  a  1:1  mapping.  Although  the  present  analytic 
scheme  has  proved  successful,  changes  in  engineering  practice  in  the 
aircraft  industry  may  make  surface  fitting  a  more  attractive  approach  in 
the  future. 

It  should  be  noted  that  as  the  SCT  program  is  presently  structured,  a 
surface  spline  would  be  very  difficult.  The.  main  computational  program 
is  expecting  to  map  body  data  in  a  cylindrical  coordinate  system  Into  the 
computational  coordinate  system.  The  surface  fitting  of  a  shape,  such  as 
the  lifting  body,  with  a  surface  spline  in  a  cylindrical  coordinate  system 
is  very  difficult  as  experience  with  the  DeJarnette  code,  which  uses  a 
surface  fit,  has  pointed  up  (Section  III-3  and  Reference  7).  As  a  general 
rule.  If  surface  splines  are  desired  for  fitting,  the  user  should  select  a 
coordinate  system  based  on  the  type  of  surface  to  be  fitted.  This,  of  course, 
can  be  done  if  the  mapping  is  not  1:1.  However,  Solomon  (Reference  6)  has 
observed  computational  difficulties  in  this  program  with  nonuniformiy 
distributed  body  points.  Points  which  are  uniformly  distributed  in  the  mapped 
space  to  meet  this  requirement ,  may  not  be  appropriately  spaced  in  the 
physical  space.  These  problems  may  be  resolved  In  the  future. 

A  program  to  develop  a  design  for  a  new  high-speed  flight  vehicle  for 
flight  test  experiments  has  resulted  in  a  large  amount  of  experimental  data 
(Reference  38)  about  one  candidate  shape,  the  X-24C-10D  (Figure  1). 
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Although  some  care  was  necessary  to  accomplish  them,  successful  flow 
fields  were  calculated  for  angles  of  attack  from  0°  to  6°  in  2°  increments 
at  a  freestream  Mach  number  of  5.95.  The  model  shape  was  approximated 
by  the  shape  shown  in  Figure  lb.  Results  of  these  studies  were  presented 
by  Neumann  at  the  AIAA  16th  Aerospace  Science  Meeting  (Reference  39) . 

Great  care  was  exercised  in  the  selection  of  body  points.  As  the 
angle  of  attack  increases,  the  low  pressure  on  the  upper  surface  causes 
the  velocity  in  the  circumferential  direction  to  increase,  resulting, 
eventually,  in  the  axial  velocity  component  becoming  subsonic  along  the 
leading  edge.  To  avoid  this  region  points  were  permitted  in  the  region 
75°  <<}><_  90°.  Solomon  (Reference  6)  states  that,  as  a  minimum,  when 
nonuniform  spacing  is  used,  the  second  differences  must  be  fairly 
uniform  and  a  preprocessor  is  supplied  to  accomplish  this  arrangement  of 
the  points.  Table  II  shows  the  comparison  of  a  desired  <j>-spacing  and 
the  ^-spacing  after  smoothing.  Because  of  this  smoothing,  a  redistribution 
of  points  with  increased  spacing  in  some  region  usually  required  a  point 
be  dropped.  It  should  also  be  noted  that,  if  strong  shocks  are  forming, 
both  the  radial  and  circumferential  spacing  must  be  increased  to  permit 
the  pressure  rise  to  be  spread  over  a  larger  area.  The  presence  of  the 
strong  shock  structure  caused  numerical  instabilities  which  forced  selection 
of  first-order  accuracy  rather  than  second-order  accuracy  for  the  body 
point  calculations.  It  might  be  suspected  that  entropy  smoothing  would 
help  in  this  case,  but  this  was  not  observed.  All  runs  for  the  X-24C-10D 
study  were  made  with  7  radial  points  informally  spaced  between  the  shock 
and  body;  the  6°  angle  of  attack  case  used  only  13  circumferential  points 
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instead  of  the  15  circumferential  points  used  at  the  lower  angles  of 
attack.  The  actual  circumferential  values  used  for  these  runs  are 
shown  in  Table  III. 

The  default  value  of  MOD  3  for  the  calculation  of  the  body  points  was 
used  (see  Section  IV-2).  The  body  description  caused  the  program  to  see 
a  longitudinal  discontinuity  which  resulted  in  the  automatic  selection 
of  MOD  0.  On  some  of  the  runs,  this  selection  resulted  in  program  failure 
some  10  to  15  nose  radii  downstream  of  the  discontinuity,  even  though  the 
highest  angle  of  attack  was  only  6°.  However,  stopping  the  calculation 
about  4  nose  radii  downstream  from  the  discontinuity  and  restarting  with 
the  MOD  3  option  again  permitted  the  solution  to  continue. 

Even  with  these  limitations,  the  comparisons  of  numerical  surface 
pressures  and  experimental  values  as  shown  in  Figures  18  and  19  are  good. 
Even  on  the  leeside  of  the  model,  where  we  would  not  expect  good 
agreement  at  angles  of  attack,  the  agreement  is  reasonable.  The  leading 
edges  show  the  least  agreement  as  might  be  expected. 

Surface  pressure  is  not  an  adequate  standard  on  which  to  base  the 
validity  of  a  flow  field  program.  Unfortunately,  pitot  pressure  surveys, 
flow  angularity  measurements,  and  shock-shape  studies  were  not  available. 
There  is  a  need  to  use  a  simpler  shape  displaying  the  basic  characteristics 
of  the  X-24C-10D  model  with  more  flow-field  probing  to  validate  the  results 
of  this  program.  Such  a  study  is  being  initiated. 

The  results  of  the  SCT  with  the  general  body  description  must  be 
examined  with  a  critical  eye.  For  example,  note  the  bow  shock  shown  in 
Figures  20  and  21.  Model  Station  (MS)  5.55  appears  quite  reasonable,  but 
notice  the  "dimple”  at  <J)  =  110°  at  MS  9.85,  This  is  probably  the  inter- 
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section  of  an  imbedded  shock  from  the  leeside  of  the  leading  edge  of  the 
model  and  the  bow  shock.  In  the  same  manner,  the  "dimple"  at  <p  =  75° 
on  MS  14.25  is  probably  an  imbedded  windward  shock  intersection.  But 
what  about  the  behavior  of  the  windward  shock  at  MS  20.52?  Is  this 
caused  by  program  instability  or  shock  structure?  Additional  experimental 
information  and  a  more  detailed  examination  of  the  computed  flow  field 
will  be  required  to  resolve  this  question. 

Although  it  is  not  easy  to  see  in  Figure  20,  there  is  an  inflection 
of  the  windward  shock  in  the  symmetry  plane.  In  Figure  22  the  slope  of 
the  shock  in  the  axial  direction  is  plotted  and  a  strong  inflection  at 
MS  5.86  can  be  seen  as  well  as  a  weaker  one  at  MS  12.51.  The  first 
inflection  had  a  significant  effect  on  the  program,  resulting  in  clearly 
observable  instabilities  in  the  flow  which  eventually  resulted  in  negative 
entropies  being  calculated,  a  clear  warning  of  danger.  No  failures  were 
observed  in  this  region,  however. 

The  first  inflection  is  easily  visible  in  Schlieren  photographs 
taken  during  the  X-24C-10D  wind  tunnel  tests  at  AEDC  (Reference  38). 

The  second  inflection  is  not  apparent. 

In  Figure  23  the  slope  of  the  shock  in  the  axial  direction  is  plotted 
as  a  function  of  circumferential  angle  for  three  model  stations.  The 
slopes  from  $  =  50°  to  <f)  =  120°  are  certainly  suspect  at  MS  9.85,  but 
the  behavior  damps  out  at  MS  20.52. 

Figure  24  shows  the  slope  of  the  shock  in  the  circumferential 
direction.  The  slopes  support  the  "dimples'1  shown  in  Figure  21,  but 
offer  few  clues  as  to  the  validity  of  the  technique. 
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The  foregoing  plots  (Figures  20-24)  are  presented  as  examples  of 
program  parameters  which  must  be  considered  in  validating  the  results 
of  the  program.  They  are  not  an  exclusive  set;  such  questions  as  Mach 
number  behavior,  entropy,  and  flow  angle  must  be  examined  also.  It 
should  also  be  noted  that  disturbances  in  the  flow  field  do  not 
necessarily  cause  problems  on  the  surface  of  the  body.  The  disturbance 
must  propagate  within  the  Mach  cone  toward  the  body,  striking  if  far 
downstream  of  the  original  disturbance.  By  this  time,  the  numerical 
differencing  may  well  have  damped  the  disturbance  to  a  point  where 
it  is  negligible. 
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SECTION  V 

HEAT  TRANSFER  CALCULATIONS  USING  SCT  PRESSURES  AND  THE  DEJARNETTE  PROGRAM 

A  major  goal  of  this  effort  has  been  to  show  numerically  generated 
surface  pressures  can  be  used  in  the  heat  transfer  code  by  DeJarnette 
(Reference  4)  to  produce  reliable  heat  transfer  rates.  We  have  found 
two  serious  limitations  to  this  approach;  there  is  a  limit  on  the 
generality  of  the  geometry  the  DeJarnette  code  will  accept  (the  X-24C-10D 
has  not  been  programmed  in  the  DeJarnette  code  because  of  this  difficulty) 
and  the  calculated  pressures  must  be  surface  fit  in  a  manner  that  permits 
the  first  and  second  derivatives  cf  the  pressure  with  respect  to  the  axial 
and  circumferential  coordinates  to  be  obtained.  These  pressures  do  not 
constitute  a  well-posed  problem  mathematically  and  the  resultant 
derivatives  are  not  necessarily  accurate.  It  is  therefore  necessary  to 
examine  the  results  of  the  surface  fit  and  the  derivatives  with  some  care 
to  insure  they  are  physically  realistic.  No  pressure  problems  were 
encountered  with  the  shape  selected  for  this  study,  an  ogive-cylinder. 

The  ogive-cylinder  used  for  this  study  is  shown  in  Figure  15.  Surface 
pressure  data  computed  by  the  SCT  Program  was  saved  and  punched  into  cards 
in  the  format  used  by  the  DeJarnette  code.  Comparisons  between  the  numerical 
and  experimental  values  were  made  for  raod-1  station  (MS)  47.5  at  angles  of 
attack  of  0°,  4°,  3°  and  12°  at  a  freestteam  Mach  number  of  5,95.  The 
wind  tunnel  model  contained  75  thermocouples  in  one  quadrant  of  the  model 
at  MS  47.5.  The  model  was  injected  into  the  tunnel,  pitched  to  the  desired 
angle  of  attack,  and  rolled  180°  in  30°  increments  with  a  set  of  heat 
transfer  data  taken  at  each  condition  for  a  total  of  525  bits  of 
experimental  data.  The  experimental  Stanton  numbers  were  calculated  for 
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each  of  these  bits  of  data  using  freestream  velocity,  density,  and  total 
temperature  rather  than  the  adiabatic  wall  temperature  and  boundary  layer 
edge  conditions  which  are  unknown. 

Pressure  data  were  also  taken  at  MS  47.5  at  18  points  around  the 
body.  Neumann  (Reference  37)  found  that  if  the  surface  pressure  was 
normalised  by  the  freestream  static  pressure  and  plotted  against  the 
tunnel-based  Stanton  number  on  a  log-log  plot,  the  result  was  a  straight 
line  of  the  form 

log  St„  »  1.081  log (p/poo)  +  0.5922.  (7) 

Figure  25  shows  the  experimentally  obtained  Stanton  numbers,  the 
Stanton  numbers  based  on  the  linear  relationship  above  and  the  Stanton 
numbers  predicted  by  the  DeJarnette  Program.  As  expected,  in  the  regions 
cf  flow  separation  on  the  leeside  of  the  model,  the  Stanton  numbers 
are  not  well  matched  by  the  DeJarnette  values.  On  the  windward  surface, 
where  good  agreement  was  expected,  the  DeJarnette  values  are  consistently 
lower  than  the  experimental  values.  It  should  be  noted  that  the 
numerical  pressures  obtained  by  the  SCT,  and  on  which  the  DeJarnette. 
values  are  based,  are  also  lower  than  experimental. 

To  insure  the  observed  differences  are  not  due  to  the  manner  in  which 
the  Stanton  number  is  calculated,  the  heat  transfer  rates  were  compared 
directly  in  Figure  26.  Only  the  4°  angle  of  attack  case  was  examined 
and  not  all  available  experimental  data  points  were  plotted,  but  it  is  clear 
the  pattern  is  the  same  as  that  observed  for  the  Stanton  numbers.  In 
Figure  27  the  experimental  and  numerical  Stanton  numbers  on  the  windward 
centerline  are  compared  as  a  function  of  angle  of  attack.  Also  shown  in 
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Figure  27  are  the  numerical  surface  pressures  and  the  experimental  surface 
pressures  for  angles  of  attack  of  4°,  8°  and  10°  (no  surface  pressure 
data  was  taken  experimentally  at  12°).  The  pressures  at  0°  angle  of  attack 
are  not  shown,  but  are  almost  identical,  0,85  for  the  numerical  case 
and  0.86  for  the  experimental.  The  differences  in  pressure  are  not 
sufficient  to  explain  the  resultant  Stanton  number  differences,  but 
they  are  in  the  correct  direction. 

A  better  geometry  and  pressure  fitting  program  must  be  developed  for 
the  DeJarnette  code  to  permit  heat  transfer  studies  on  more  complex 
shapes  such  as  the  X-24C-10D.  A  program  with  McDonnell  Douglas 
Astronautics,  East  has  been  started  to  meet  this  requirement. 
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SECTION  VI 
CONCLUSIONS 

The  original  objective  of  this  study  was  to  demonstrate  that  the 
heat  transfer  portion  of  the  HABP  is  useful  for  lifting  reentry  shapes. 
This  objective  was  realized  only  in  a  limited  manner,  as  the  restrictions 
imposed  on  streamline  tracking  limit  the  usefulness  of  the  technique. 

Heat  transfer  rates  obtained  with  this  method  for  the  flat  portion 
of  the  X-24C-10D  side  panel  compare  well  with  data. 

Within  the  restrictions  and  limitations  of  the  DeJarnette  program 
it  represents  a  fairly  simple,  relatively  accurate  computer  program  to 
calculate  aerodynamic  heating  rates  over  fairly  general  three-dimensional 
bodies  and  is  available  to  AFFDL .  The  program  is  fast,  requiring  on 
the  order  of  seven  seconds  to  complete  a  screamline  for  typical  cases 
described  in  this  report.  In  general,  it  requires  that  the  surface 
pressure  distribution  be  obtained  independently;  however,  it  has  the 
capability  of  computing  a  modified  Newtonian  pressure  distribution. 

An  improved  description  of  the  inviscid  flow  field  was  required  to 
obtain  the  improved  pressure  distribution.  The  NSWC  Shock- Capturing 
Technique  was  obtained  and  made  operational  to  meet  this  requirement. 

The  code  was  used  to  calculate  the  inviscid  flow  about  biconics,  ogive- 
cylinders  and  a  lifting  body.  Because  of  limits  in  the  existing 
DeJarnette  code,  it  was  not  possible  to  use  the  DeJarnette  code  with 
the  lifting  body;  however,  pressures  calculated  by  the  Inviscid  program 
for  the  ogive-cylinder  were  put  in  the  DeJarnette  code.  The  resultant 
heat  transfer  rates  compared  well  with  experiment. 
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In  summary,  techniques  are  available  to  calculate  the  surface 
pressure  distribution,  streamline  pattern  and  heat  transfer  rate  to 
a  large  class  of  axisymmetric  bodies  at  angle  of  attack.  More  limited 
capabilities  have  been  demonstrated  for  bodies  with  a  single  plane  of 
symmetry  at  angle  of  attack. 
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APPENDIX  A 

BODY  PROCRAM  FOR  OGIVE-CYLINDER 

In  the  three  dimensional  shock  capturing  technique,  the  body  radius 
and  derivatives  required  by  the  main  program  are  calculated  in  a  subroutine 
called  BODY  and  transmitted  back  to  the  calling  program  through  a  common 
block.  The  values  of  (f>-angle  for  which  calculations  are  required  and 
the  value  of  the  longitudinal  station  are  also  transmitted  through  common. 
Thus,  only  one  calling  argument,  the  index  number  of  the  cross  sectional 
plane  being  calculated,  is  needed.  The  original  BODY  subroutine  can 
be  replaced  by  any  new  subroutine  meeting  these  requirements. 

An  ogive-cylinder  may  be  described  by  giving  a  radius  for  the 
ogive  portion  or  the  length  of  the  ogive  portion  plus  the  radius  of  the 
cylinder  as  shown  in  Figure  A-l.  The  BODY  program  must  calculate  the 
local  body  radius  normal  to  the  longitudinal  axis  along  with  the  first 
and  second  derivatives  with  respect  to  circumferential  angle  and 
longitudinal  distance  plus  the  cross  derivative.  The  calculated  shape 
must  have  a  blunt  nose  and  all  lengths  must  be  normalized  by  the  nose  radius. 

The  subroutine  written  for  this  purpose  expects  the  length  of  the 
ogive  portion  and  the  radius  of  the  cylinder  portion,  both  normalized  by 
the  nose  radius,  to  be  given  as  input.  Then,  by  Figure  A-l,  the  radius  of 


the  ogive  is  given  by 


R  -  0.5(S  2H  2)/R 
og  og  cy  cy 


(A-l) 


and  the  slope  at  the  sphere-ogive  junction  is  given  by 

a  =  acos  (k  -  R  )/(R  -  1)  (A-2) 

og  cy  og 
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and  Che  length  of  the  blunted  ogive  is  given  by 

L  ■  (R  -  1.)  sin  a  +  1.  (A-3) 

ag  og 

At  any  value  of  z  the  body  radius  is  given  by 


®B  -  : 

where 

R. 

temp 
D 

temp 

The  derivatives  are  given  by 

dRg/d<f)  -  d^/dtj^  -  d^Rg/ dzdtj)  -  0  (A-6) 

dR_ ! dz  =  D  /R 
t  temp  temp 

dV/dz2  -  -(R  +  D  )  (dR  /  dz)  /R,  2  (A-7) 

n  temp  temp  B  temp 

These  values  are  calculated  once  for  each  body  station  when  z  is  less 

than  L  .  Subsequent  calls  to  BODY  result  in  an  immediate  return  to  the 
og 

calling  program.  When  z  is  greater  than  L^,  the  cylinder  radius  is 
returned  with  all  derivatives  set  equal  to  zero. 


1  -  (R  -  R  ) 

cy  og  temp 


(A-4) 


4 


2  2 
R  -  D„ 
og  temp 


(A-5) 


■  L  -  z 
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APPENDIX  B 

BODY  PROGRAM  FOR  A  GENERAL  BODY 

In  the  three-dimensional  shock-capturing  technique  the  body  radius 
and  derivatives  required  by  the  main  program  are  calculated  in  a  sub¬ 
routine  called  BODY  and  transmitted  back  to  the  calling  program  through 
a  common  block.  The  values  of  ({'■■angle  for  which  calculations  are 
required  and  the  value  of  the  longitudinal  station  are  also  transmitted 
through  common.  Thus,  only  one  calling  argument,  the  index  number  of 
the  cross-sectional  plane  being  calculated,  is  needed.  The  original 
BODY  subroutine  can  by  replaced  by  any  new  subroutine  meeting  these 
requirements . 

The  initial  shape  considered  for  the  new  body  routine  is  shown  in 
Figure  B-l.  The  body  may  also  have  flat  surfaces  on  the  top  and  bottom. 
While  the  present  version  of  this  BODY  program  is  more  general  than  the 
original  routine  and  the  ogive-cylinder  routine,  there  are  some  specific 
restrictions . 

The  shape  may  be  divided  into  any  number  of  longitudinal  sections, 
but  each  section  may  have  no  more  than  five  segments  in  the  circumferential 
direction.  These  segments  may  be  arcs  of  circles  or  straight  lines  and 
may  have  longitudinal  dependence. 

Prior  to  starting  each  new  segment  of  the  body,  the  program  computes 
the  maximum  (ji-augle  for  that  segment.  Simple  geometry  is  used  to  locate 
this  angle;  for  example,  if  the  current  segment  and  the  next  segment  are 
both  arcs  of  circles,  then  the  current  segment  will  end  when  its  radius 
passes  through  the  center  of  the  next  arc.  Going  from  a  circular  arc  to 
a  straight  line,  the  arc  ends  when  the  slope  of  the  arc  equals  the  slope 
of  the  straight  segment.  As  can  be  seen,  this  and  similar  procedures  used 
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for  the  other  possible  cases  greatly  simplifies  the  body  description. 
Because  not  ail  bodies  have  continuous  slopes,  it  may  be  wise  to 
supply  the  maximum  <}>  angle  for  each  segment. 

A  circle  is  described  by  giving  the  coefficients  to  specify  the 
locus  of  the  center  and  the  radius,  thus 

x, -AZ2  +  BZ+C 

ax  xx 

y.  «AZ2  +  BZ+C  (B— 1) 

k  y  y  y 

R  =AZ2+BZ+C 
or  r  r 

A  straight  line  intersecting  the  plane  of  symmetry  is  uniquely  described 
by  the  slope  relative  to  the  horizontal  axis  and  its  intercept  point  on 
the  vertical  axis 

x=AZ2+BZ+C 
n  x  xx 

A  general  straight  line,  not  intersecting  the  plane  of  symmetry,  is 
uniquely  defined  by  its  adjacent  segments  assuming  slopes  are  continuous 
about  the  body.  However,  to  simplify  calculations,  the  program  expects 
the  slope  to  he  given. 

The  first  segment  at  each  station  may  be  either  a  straight  line  or 
a  circular  arc.  If  the  first  segment  is  a  straight  line,  the  second 
segment  must  be  a  general  circular  arc.  The  body  must  have  a  plane  of 
symmetry  and  no  logic  is  included  to  handle  angles  greater  than  180°. 

Having  decided  to  use  an  analytic  approach  to  the  body  geometry,  the 
problem  to  be  solved  is:  Given  a  body  described  in  a  cartesian  system  and 
given  an  angle  <J>  in  the  circumferential  direction,  find  the  body  radius 
and  its  first  and  second  derivatives  with  respect  to  Z  and  <j>  plus  the 
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cross  derivative.  The  body  radius  can  be  quickly  calculated  by  the  Law 
of  Sines  using  Figure  B-2  for  a  circular  arc  crossing  the  plane  of  symmetry. 

R^  «  Rq  sin  6/sin  ij > 


where 


0  -  180  -  ip  -  Y 
Y  “  asin  (Bgsinp/Ro) 

Rq  =  R(Z),  Bg  -  B(Z),  which  makes  Y=Y(Z),  ip  =  ip(p),  and  R^  “  f(p,Z). 

The  derivatives  of  with  respect  to  Z  arid  <j>  are  straight  forward,  but 

long.  They  are  listed  here  with  some  intermediate  derivatives  left 
unstated. 


dRB  d0 

_  _  R  sin  i)  cos  0  _ 4-  cos  ip  sin  G 

d<j)  °  dp 

sin2p 


(B-3) 


2 

— -j—  Ro  f  sin2ip  (-simp  sin9  (~)2  +  simp  c.os@  ^-4  (B-4) 

d<P  JTnH  L  d<J> 

+  simp  sinG)  -(simp  cosG  ^4  +  cosip  sin9)  (-2sinipcosip)j 


dZ 


*■  sin0  dR  R  cos6  dY 
o  —  o 


(B-5) 


dZ 


dZ 


dZ2 


sin  i p 


,  2 

d  Ro  sin  G-2  dR°  cos  0  dY  -  R  sin  0 
dZ7™  dZ  dZ  ° 


(S)2 


(B-6) 


A,  \ 


-  2R  cos  0  d  Y  |  /  sin  ip 
0  dZ 7 


d  -  (cos  6  d6_  dRn  +  R  sin  0  d0_  dY  -  R  cos  0  d2Y  V  <B-7) 

dp  dZ  dZd 


dZd 


P 


dip 


dZ 


sin  ip  +  (sin  0  dR  -  R  cos  0  dY  \cos  \{i/sin  ^ 

dT  °  dz 
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Note  all  derivatives  become  undefined  at  4*  =  0°  an^  180°.  A  check  of  the 
numerators  shows  these  forms  are  all  indefinite  (0/0)  and  by  L'Hospital's 
Rule  we  find 


%  " 

R  -  B 
o  s 

(B-8) 

-  d2RB  = 

d\  -0 

(B-9) 

dc}> 

dtj>2 

dZd(p 

dRfi 

1 

P^0j 

11 

dB 

s 

(B-10) 

dZ 

dZ 

dZ 

2 

,2 

dRB 

=  d2R  • 

dB 

s 

(B-ll) 

diT 

o 

d&r 

dZ2' 

When  calculating  for  a  general  circular  arc  as  shown  in  Figure  B-3,  we 

let  B  be  the  distance  to  the  center  of  the  circular  arc  and  <j>  be  the 
s  m 

angle  of  B  relative  to  the  axis  of  symmetry,  then  ij>  can  be  defined  as 
s 

shown.  The  derivatives  for  <f>  are  the  same  as  above,  but  the  Z-derivatives 
are  now  given  by 

dRB  "  sin  *  <sin  0  +  R0  cos  6  d6j  -  Rq  sin  0CQa  y  dc(^  (B-12) 

dZ  dZ  dZ  dZ~ 

sin2^ 


where 


da 

dZ 


dR„  cos  ^  dl^m 

dT  'd“ 


(B-13) 


tip  “  R  cos  0  d0  +  sin  0  d  R  -  R  (sin  e/d(A  " 

dZ  °  dZ~  dZ^  °  VdzV 

cos  6  d20  ^  (B-14) 

dZ 7/ 

2 

jo  g  jo  ^  dQ 

da  “  dR  cos  m  -  R_(sin  mi  -  cos  8  m\  (B-15) 

dz  —■  d r  VdzV  “3z 

aZ 
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Again  the  question  of  limits  as  <j)  -*■  0°  must  be  addressed.  The  first 
derivative  with  respect  to  remains  equal  to  zero,  but  the  second 
derivative  is  an  analytic  singularity  not  equal  to  zero  and,  at  the  present 
time,  not  properly  defined  in  the  program.  The  Z-derivatives  limits 
are  now  given  by 


dR-  dB  dR 

o  m  s  +  0 

(B-16) 

dZ  dZ  dZ 

2  2  2 
d  L,  d  B  d  R 

H  _  S  j_  O 

(3-17) 

dz2“  ~  dz2“  +  dZ? 

To  calculate  the  required  values  on  a  straight  bottom  (a  straight  top 
uses  the  same  logic),  the  construction  shown  in  Figure  B-4  is  used.  Note 
that  the  length  along  the  x-axis  is  given  and  the  slope  of  the  bottom 
relative  to  tne  horizontal  is  known.  The  radius  and  derivatives  are 
given  by 

Rg  ■  Bg  sin  0/sin  if 


6  “  ir/ 2-slope 

ip  “  ir/2  +  slope  -  cj> 


2 

dRp  “  B^  sin©  cosif  d8  /  sin  if 


d2RB  -  sine  d2Bg 

dz*-  'sin*  dZ^ 

2 

d  Rg  “  sin0  cosif  dBg  dB^ 

-11  dtt> 

sin'ty 


(B-18) 

(B-19) 

(B-20) 

(B-21) 

(B-22) 
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The  last  general  case  Is  a  general  straight  line  whose  slope  is 
measured  from  the  vertical  as  shown  in  Figure  B-5.  For  this  case 
is  a  known  length  and  makes  an  angle  $  relative  to  the  plane  of 
symmetry.  Then  G  is  ^  -  Slope  and 

B  =*  B  sin  0  (B-24) 

ms 

ip  =  n/2  +  slope  -  <j> 

From  these  values  we  have 


h  =  V003* 

(B-25) 

dRB  =  -  Bm  sin  ^ 

(B-26) 

d<(>  cos2i|i 

d2R_  *>  B 

— (i  + 

d< p*  cost|i 

2  tan2  i|j) 

(B-27) 

dVi  dBm 
dT  cos*  az~ 

(B-28) 

dV  -  1  d2B 

B  — — ,  m 

dZT  ^  dz?~ 

(B-29) 

2 

d  It,  »  -  sin  ip 

dZd<|) 

dB 

m 

dZ 

(B-30) 

In  the  case  of  the  ogive-cylinder,  all  body  radii  and  derivatives 
were  easily  obtained  to  check  the  computer  code.  In  this  case,  the 
derivatives  were  not  easily  computed.  To  check  the  values  computed  by 
the  new  body  code  a  test  routine  was  used  to  calculate  the  required 
values  at  2°  intervals  around  the  body  and  any  desired  number  of 
body  stations.  For  each  section  of  the  body,  four  stations  0.5 
nose  radii  apart  were  computed.  A  small  number  of  body  radii  were 
hand-calculated,  then  these  and  the  radii  computed  by  the  test 
program  were  carefully  plotted. 
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Once  the  body  radii  were  validated,  they  were  used  in  a  four-point 
numerical  differentiating  formula  to  calculate  the  <{>  and  Z  derivatives. 
The  differentiating  formula  used  was  (Reference  45) 

dB1/dK  ■=  (1/h)  (-2Bq  -  3BX  +  6B2  -  B3)  (B-31) 

where  K  was  <f>  or  Z  and  h  the  interval  between  B^  and  End  points 

were  checked  with 

dB0/dK  =  (1/h)  (~11B0  +  18B1  -  9B2  +  2B3)  (B-32) 

and 

dB3/dK  =  (1/h)  (-2B0  +  9B1  -  18B 2  +  11B3) .  (B-33) 

When  the  first  derivatives  had  been  checked,  the  same  pattern  was 
repeated  to  check  the  second  derivatives  and  cross  derivative.  It  is 
recommended  that  each  new  shape  be  tested  in  this  manner  prior  to  use 
in  the  SCT. 
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TABLE  II 

COMPARISON  OF  DESIRED  AND  SMOOTHED  PHI  ANGLES 


Desired 

Smoothed 

Desired 

Smoothed 

0.0 

0. 

115.0 

111.05 

23.0 

22.97 

135.0 

130.46 

46.0 

45.31 

155.0 

148.24 

69.0 

68.33 

170.0 

164.56 

93.0 

90.21 

180.0 

180.0 

TABLE 

;  in 

LOCATION  OF  CIRCUMFERENTIAL 

PLANES  FOR  X- 

24C-10D  CASES 

Values  of  <p 

for  Angles 

of  Attack  of 

0°,  2°,  4° 

0.00 

30.44 

64.85 

107.63 

151.89 

10.04 

41.14 

78.27 

122.70 

166.01 

20.14 

52.52 

92.66 

137.49 

180.00 

Values 

of  Angles  of  Attack 

.  of  6° 

0.00 

42.42 

90.07 

136.66 

180.00 

13.61 

57.91 

106.93 

151.33 

27.66 

73.92 

121.58 

165.73 

! 
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Figure  4.  Plot  of  Streamline  Angle  Vs  Angle  of  Attack 


Orthogonal  Coordinate  System  and  Body  Orientation  for  the 
DeOarnetta  Program 
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Figure  7.  25°  Cone  Heat-Transfer  Distribution 
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Figure  10.  Continued 
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Figure  12.  00°  Delta  Wing  Spanwise  Heat-Transfer  Distribution 
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Figure  13.  10°  Cone  Windward  Streamline  Heat-Transfer  Distribution 
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Figure  16.  Surface  Pressures  on  Ogive-Cylinder  at  Model  S  : 
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Figure  18.  Continued 
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Figure  19.  Surface  Pressures  on  the  X24C-10D  at  6°  Angle  of  Attack 
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at  Three  Body  Stations 
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Figure  24.  Concluded 
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Figure  25.  Comparison  of  Hunerica''  and  Experimental  Station  Number: 
on  the  Ogive-Cylinder 
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Numerical  and  Experimental  Heat  Transfer  to  the 
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Comparison  of  Numerical  arid  Experimental  Pressures  and 
Stanton  Numbers  on  the  Ogive-Cylinder  as  a  Function  of 
Angle  of  Attack 
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Figure  B-5.  Construction  Used  for  a  General  Straight  Segment 
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